Équation de propagation de Bessel, construction des solutions 


L'équation de propagation de Bessel ou en anglais < Bessel Wave Equation > est paradoxalement 
une fonction peu étudiée, en réalité à part l'ouvrage de P.Moon et D.E.Spencer « Field Theory 
Handbook », il en est peu fait mention dans des articles de recherche, ou alors de façon plus ou 
moins directe, notamment qui établissent en préambule que ces fonctions peuvent se construire 
formellement à l'aide des fonctions de Whittaker M et W, (voir les trois références suivantes : 
H.Hochstadt « The Functions of Mathematical Physics» chapitre 7, 1957 H.Hochstadt 
"ADDITION THEOREMS FOR SOLUTIONS OF THE WAVE EQUATION IN PARABOLIC 
COORDINATES" Pacific Journal of Mathematics Vol 7, n*3, A.Erdelyi et H.Bateman 
« HIGHER TRANSCENDENTAL FUNCTIONS VOLUME П, Chapitre ҮШ, Functions of the 
Parabolic Cylinder and of the Paraboloid of Revolution»). D'autres références comme dans le livre 
des auteurs M.Willatzen, L.C.Lew.Yan.Voon < Separable Boundary-Value Problems in Physics > se 
contentent de reprendre la construction de Frobënius en plagiant « P.Moon et D.E.Spencer », mais 
uniquement sur la fonction de première espèce. 


Pourtant cette équation intervient par exemple dans la séparation des variables de l'équation de 
Helmholtz dans au moins deux systèmes de coordonnées : 

- parabolique cylindrique, si la solution ne dépend pas de z, comme l'équation de la chaleur sur un 
cylindre de section parabolique et de hauteur infinie 

- parabolique de révolution (ou paraboloïdal de révolution), si la solution ne dépend pas de l'angle 
de révolution, comme l'équation de la chaleur sur un paraboloïde de révolution 


Certes il convient de se méfier de certaines approches de P.Moon et D.E.Spencer dans « Field 
Theory Handbook >, non qu'elles soient fausses, bien au contraire mais qui laisse apparaítre de 
nouvelles dénominations de fonctions alors même que ces fonctions avaient à l'époque du livre déjà 
été amplement étudiées. Je n'en veux pour preuve que la fameuse équation de propagation de 
Legendre, donnant possiblement des fonctions d'onde de Legendre alors même que l'équation se 
nomme plus communément équation d'ondes sphéroïdales donnant lieu aux fonctions d'ondes 
sphéroïdales, très largement étudiées par Meixner et Schafke dans « Mathieu Functions and 
Spheroidal Functions and Their Mathematical Foundations ». Toujours est-il qu'après des 
recherches de plusieurs mois dans la littérature scientifique, force est de constater de façon 
surprenante que les fonctions d'onde de Bessel sont en quelque sorte une « terra incognita », sauf à 
constater que l'équation est proche de l'équation de Whittaker, et qu'elle peut servir à la 
construction formelle de solutions. Il n'en demeure pas moins que l'analogie réalisée par P.Moon et 
D.E.Spencer avec les fonctions de Bessel est probablement féconde car elle nous conduit à tenter 
de construire les solutions « à la manière » des fonctions de Bessel. 


Le plan d'étude est assez classique par la construction des solutions indépendantes par la méthode 
de Frobënius, puis les relations des solutions avec les fonctions de Whittaker M et W, et également 
avec les fonctions d'onde de Coulomb. Il convient d'établir les diverses formes de fonctions d'onde 
de Bessel à l'image des fonctions usuelles de Bessel J et Y, et de Hankel H(1) et H(2). Il convient 
également de donner autant que faire se peut, les solutions de l'équation d'onde de Bessel modifiée 
et solutions correspondantes I et K. 


L'équation de propagation de Bessel prend la forme suivante : 
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Et l'on pose par convention solution Z(x)= AJ, (a. В,х)+ B,J. (e. B.x) 


On lui donne le nom de Bessel car elle en constitue une généralisation et redonne l'équation usuelle 
de Bessel en posant a=0. 


L'équation de propagation de Bessel modifiée prend la forme suivante : 
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Séparation des variables et équation de propagation de Bessel en coordonnées paraboloïdales de 
révolution 
Comme nous le disions c'est dans le système de coordonnées paraboloïdales de révolution (ou 


paraboliques de révolution) que l'on rencontre l'équation de propagation de Bessel. Dans ce 
système le Laplacien prend la forme suivante : 
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Par séparation de l'équation de Helmholtz, il vient : 
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Avec les autres variables paraboliques, nous avons : 


n-42€ vog = х= uv Cos(o) »-uvSin(p) z= =" 


2 
(nd ag reg fn) 4r iP) B (au) 
u du du u 
ct LO [at + B? EE 02 f,(v)= AJ, (a, B,v)+ B.J | (o. B.v) 
v dv dv v 
` Tri, viel 
9 


Les deux premiéres équations différentielles sont bien de la catégorie des équations de propagation 
de Bessel. Lorsque l'on envisage un probléme de diffusion, les équations sont de la méme forme. 


Séparation des variables en coordonnées paraboliques cylindrique sans dépendance à la hauteur z 


Le systéme de coordonnées paraboliques polaire (désignant le systéme de coordonnées restreint au 


plan z-0: „ _ u-v? Е . Ге Laplacien en coordonnées parabolique 2D (u,v) s'écrit 
х=а S y=auv 
2 2 
sous la forme ` AT(u,v,t) = — | : Ô сЗа д 0 
a (u +v ) ди ду 

2 2 

L'équation de diffusion s'écrit donc : l 0 T(u,v,t) + OT(u,v;t) | 1 ôT(u,v,t) _ 0 
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Par séparation des variables la solution de l'équation s'écrit : 
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Si l'on poursuit la séparation des variables sur les deux coordonnées : 
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Ce sont des équations assimilables à l'équation de propagation des ondes de Bessel 
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Lien avec les fonctions de Whittaker 


Nous allons voir qu'en transformant l'équation des ondes de Bessel, on retombe sur une équation 
plus connu, l'équation de Whittaker. 


Partons de l'équation de propagation de Bessel, pour changer de variable et de fonction 
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Cette dernière équation est à très peu près l'équation de Whittaker, donnant lieu aux solutions qui 
sont des fonctions de Whittaker de première et deuxième espèce . Un dernier changement de 
variable est réalisé : 


" a? p 1- p? 
y Où! 1 od P prono 


sr ] EL. EE E EISE ЗА Zu: 
4 4а 4 


Les solutions de cette dernière équation se trouvent être des fonctions de Whittaker. Ce sont les 
fonctions M et W de Whittaker. 


V (v) = WhittakerM,., (v) 
V, (v) = WhittakerW, „ (v) 
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En définitive on peut établir un lien formelle entre des combinaisons linéaires des fonctions d'onde 
de Bessel et des combinaison linéaire des fonctions de Whittaker transformées comme suit : 
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L'inconvénient d'utiliser ces fonctions de Whittaker est qu'en général elles sont à valeur imaginaire 
du fait que l'argument et les parametres injectés le sont eux-mémes. Mais par bonheur on va voir 
qu'une renormalisation de la solution adéquate permet déjà pour les fonctions de première espèce 
de retrouver les fonctions d'onde de Bessel tel que définies à partir des développement de 
Frobenius. Il en est de même pour les fonctions d'onde de Bessel de deuxième espèce mais d'ordre 
non entier. 


L'analyse des liens entre les fonctions de propagation de Bessel de première espèce et les fonctions 
de Whittaker M est ici reporté de manière empirique. Mais il y a de fortes chance qu'il soit juste. 
L'approche consiste à renormaliser l'expression de Whittaker M de telle sorte que lorsque le 
paramètre a-70, on puisse retrouver les fonctions de Bessel. Pour cela on établit le développement 
limité des deux expressions autour de x=0, au premier ordre non nul : 
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Si de plus on trace le graphe de la fonction 


Lt 

x(ia)2 

partie imaginaire est toujours nulle. Autrement dit on peut former le résultat suivant, la fonction 
d'onde de Bessel de première espèce est construite à partir des fonctions de Whittaker M comme 
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Ce lien établit empiriquement, est bien confirmé par la formule de continuation analytique des 
fonctions de Whittaker M (voir NIST): 
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La propriété de continuation analytique est en fait l'application des transformations de Kummer sur 
les fonctions hypergéométriques confluentes M et U, soit : 
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Transformation de Kummer 
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Donc les deux fonctions ont même valeur pour x=0 et de plus sont toutes deux solutions de 
l'équation différentielle. Par la suite nous verrons également que ces deux fonctions on exactement 
le même comportement asymptotique. Elles sont donc égales et il vient : 
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Pour peu que le paramètre soit un réel, l'expression est tout aussi valable avec des valeurs 
négatives de l'ordre. Ainsi : 
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De plus ce résultat n'est au moins jamais démentie en traçant les graphes comparatif des fonctions 
de propagation de Bessel, construite par le développement de Frobenius, et l'expression analytique 
précédente. Le décrochement des solutions de Frobenius vient de la convergence trop lente des 
séries et du besoin d'un nombre de plus en plus important de termes dans les grandes valeurs de 
l'argument : 


a=12.5; 
B=1; 
p = 1; 


—— BesselJWaveRecurrentSeriesFunction(a, B, p, x) 
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а= 1.5; 

В = 1; 

p = -0.55; 


— BesselJWaveRecurrentSeriesFunction(a, 8, p, x) 
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— BesselJWaveRecurrentSeriesFunction(a, 8, р, x) 
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Ге cas des solutions de deuxième espèce d'ordre entier semble plus compliqué (voir le calcul plus 
loin à l'aide des fonctions de Coulomb). L'expression finale est remarquablement compacte, elle 
porte sur une fonction de deuxième espèce Y légèrement modifiée que nous nommerons Y tilde : 
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La fonction Y tilde n'est pas exactement celle construite par le développement de Frobenius. Elle 
diffère d'un terme proportionnel à la fonction d'onde de Bessel de première espèce, qui s'annule 
lorsque a=0 (cas des fonctions de Bessel ordinaire). Nous allons par ailleurs chercher à déterminer 
ce terme de proportion, à l'aide d'un développement asymptotique (voir plus loin). 


Et voici quelques graphes de cette fonction, représentée avec la fonction Y de Neumann classique : 
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— (х) 
— BesselYWaveFromWhittakerWCompact(1.5`, 1, 0, x) 
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— Y3(x) 
— BesselYWaveFromWhittakerWCompact(3.5 , 1, 3, x) 
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— Y1(2x) 
— BesselYWaveFromWhittakerWCompact(1.5`, 2, 1, x) 


—— BesselYWaveFromWhittakerWCompact(3.5`, 2, 2, x) 
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Développement limité des fonctions d'onde de Bessel J d'après les fonctions de Whittaker 


Par curiosité, on peut étudier le développement limité autour de x=0 de l'expression des fonctions 
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Il est parfaitement clair qu'en tout point le développement limité coincide avec celui du 
développement de Frobenius. Nous avons pu le tester notamment jusqu'à l'ordre 20 à l'aide de 
Mathematica, ce qui nous a paru suffisant comme test. On pousse le développement limité jusqu'à 
l'ordre 18 pour la beauté du geste : 
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Fonction d'onde de Bessel modifiée de première espèce IÍ 


Par construction nous avons vu qu'avec le changement de variable x->ix, on obtenait la définition 
des fonctions d'onde de Bessel modifiées, soit : 
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fonctions de Bessel usuelles, les formules de connexion sont les suivantes (voir NIST Handbook of 
Mathematical function, Bessel Functions , formules 10.27.6 : 
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Par extension, on trouve l'expression des fonctions d'onde de Bessel modifiée L en ayant soin 
d'appliquer la deuxième transformation I, (o. В,х)= i" J „ба, B,-ix) cohérente avec la définition des 


fonctions de Whittaker M : 
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L'étude préliminaire des graphes des fonctions d'onde de Bessel modifiée de première espèce 
montre que contrairement à la fonction Í de Bessel la fonction oscille et se révèle décroissante à 
l'infini. On reproduit des graphes présentant à diverses valeur de a, ces oscillations : 
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— BessellWaveFromWhittakerM(0.1, B, р, x) 
— BessellWaveFromWhittakerM(0.2, 8, p, x) 
—— BessellWaveFromwhittakerM(0.3, B, p, x) 
— BessellWaveFromwhittakerM(0.4, 8, p, x) 
— BessellWaveFromwhittakerM(0.5, 8, p, x) 
— BessellWaveFromwhittakerM(0.6, 8, p, x) 
— BessellWaveFromwhittakerM(0.7, 8, p, x) 
— BessellWaveFromWhittakerM(0.8, 8, p, x) 
— BessellWaveFromwhittakerM(0.9, 8, p, x) 
— BessellWaveFromWhittakerM(1., B, p, x) 
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В = 0.01 
р= 2 

— BessellWaveFromWhittakerM(0.1, B, р, x) 
— BessellWaveFromWhittakerM(0.01, 8, р, х) 
— BessellWaveFromWhittakerM(0.001, 8, р, х) 


L'étude du développement limité de l'expression en fonction de 


développement limité de Frobënius, à savoir jusqu'à l'ordre 18 : 
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Lien avec les fonctions de Coulomb 


On a vu que l'équation de propagation de Bessel est équivalente à l'équation suivante : 


2 
Equation de propagation de Bessel 7"(х)+ 17 ()+ [ar +В? – E о) = 0 
х х 


2 2 2 
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B? => yv) [ se o = 0 € Equation d'onde de Coulomb 
v v 
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Donc les solutions de l'équation de propagation de Bessel sont proportionnelles aux solutions de 
l'équation d'onde de Bessel moyennant un changement de variable : 


V" (v)+ [ s H “Л = 0 < Equation d'onde de Coulomb 
v v 


| B? ax B? a x? 
Sol АОСУ СЕ оао EE 
olution ` (n v) TU, v)2 (x) a às dÉ sl 4a 2 


La solution régulière de l'équation de propagation de Coulomb s'écrit à l'aide de fonction de 
Whittaker : 


Fin, x)= cn) iY" Wninakertd in, 1+— TE SH d'rue in, bi 1-28) 
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Lorsque l est un entier, on peut calculer le coefficient Cin), explicitement, en utilisant les propriétés 
des fonctions de I; à savoir : 
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Par ailleurs on démontre que les deux formes de la fonction régulière de Coulomb sont 
équivalentes à l'aide des propriétés de la transformation de Kummer sur la fonction 
hypergéométrique confluente M : 


z Ï 
= Zu 1 
WhittakerM (y, u, z) =e 222 ` Content} +u—7,1+ 22) 
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Transformation de Kummer 


ConfluentM(a,b,z)= e*ConfluentM(b-a,b,-z) 
ConfluentU(a,b,z)= z ^ ConfluentU (a -b+1,2-b,-z) 
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La deuxième solution et également la première solution de l'équation de Coulomb se construisent 
comme suit : 


ZN ig (r 
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Lien des fonctions d'onde de Coulomb de première espèce avec les fonctions d'onde de Bessel de 
première espèce 


Donnons le lien entre les solutions de première espèce de ces deux équations de propagation. 
Comme on a établit : 


2 
Ba ven! L uas] 
ia 
s жт 27 
a ? i? T(1+p) 
En injectant la solution régulière de l'équation de l'équation de Coulomb comme suit, on peut 
établit le lien formel avec la solution régulière de l'équation de propagation de Bessel : 
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Comportement asymptotique des fonctions d'onde de Coulomb pour n largement négati 


Pour déterminer le comportement asymptotique des fonctions d'onde de Coulomb en relation avec 
les fonctions d'onde de Bessel, il faut combiner deux développements, celui en x tout d'abord ( voir 
NIST Hanbook of Mathematical Functions : Coulomb Functions, formule : 33.6.1) : 


Fn.x)- Cn) Z Am)" Ап) = Aln)= 0) = Polynome en n de degré k 
k=0 
Puis les comportements asymptotiques connus des fonctions de Coulomb pour une valeur 
largement négative du paramètre n, ce qui correspond à une valeur a tendant vers 0 pour les 


fonctions d'onde de Bessel (voir NIST Hanbook of Mathematical Functions : Coulomb Functions, 
formule : 33.10.7 et 33.10.8) : 


Fa 10 DED тх), жтт du! 
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Or ces deux comportements asymptotiques des fonctions Е et С пе conviennent pas tout à fait car il 
sont définis pour une valeur x constante. Or dans notre cas en réalité c'est le produit n*x qui est 
constant. Le premier développement de la fonction régulière correspond en tout point au premier 
terme du développement en x, en effet : 
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21+2)\ 2 21+2 


> F,(n.x)= ach ZE) (2л x), zack A, Canxk(-20x] ^ 


= EXD Cos = egy 
(- 20) 
On peut donc conclure de cette remarque que le terme suivant est obligatoirement en puissance 1+2 
de x. Hélas n'ayant pas trouver un développement en x équivalent pour la fonction irréguliére G je 
ne peut réaliser le méme raisonnement. Toutefois on peut supposer qu'il en est de méme. Et il me 
semble plus exact de supposer les comportements asymptotiques comme suit : 
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Dans ces conditions nous sommes plus confortable pour appliquer les substitutions suivantes : 
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On retrouve le passage à la limite de la fonction de propagation de Bessel de première espèce 
lorsque le paramètre a->0 qui redonne une fonctions de Bessel de première espèce. 


On peut s'assurer de hypothèse sur des graphes successifs des deux expressions, comparant le 
développement de Frobenius et la fonction de Coulomb. 


Pour a=1.5, В=1, р=1 
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— BesselJWaveRecurrentSeriesFunction(a, B, p, x) 
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Pour a=4, В=1, р=3 
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— BesselJWaveRecurrentSeriesFunction(a, 8, p, x) 
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Pour а=1.5, B=1, p=-0.55 
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— BesselJWaveRecurrentSeriesFunction(a, : р, X) 
(E [trees], 


— BesselJWaveRecurrentSeriesFunction(a, 8, p, x) 
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Relation de récurrence sur les fonctions de Coulomb F et G 


Sur la fonction d'onde de Coulomb régulière F : 
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Sur Іа fonction d'onde de Coulomb G, les formules sont identiques : 
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Ces relations de récurrence appliquées aux fonctions d'onde de Bessel donne la formule suivante : 
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Le cas limite a=0, redonne normalement les fonctions de Bessel, il vient alors pour la récurrence : 


2(р-2)2(р-1\(р-3)- 8^» V, (B y)-(p-1)B*3*J, (B у) 
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En combinant les diverses relations de récurrence connues sur les fonctions de Bessel et en 
substituant les termes p-1 et p-3,il vient la méme relation : 
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Hélas on ne peut rejouer en sens inverse ce type de relation de récurrence pour en déduire pour les 
fonctions d'onde de Bessel une relation de récurrence entre plus proches voisins. 


La méthode de passage à la limite peut nous permettre une hypothèse quand à la correspondance 
de la seconde solution de l'équation de propagation de Coulomb avec la seconde solution de 
l'équation de propagation de Bessel, sachant que le passage à la limite a->0 est valide pour la 


seconde solution Ү,(о,В,х)-> Y,(Bx) : 
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Toutefois d'après le tracé des graphes, les fonctions d'onde de Bessel de deuxième espèce construite 
par le développement de Frobenius ne coincident pas tout à fait, c'est la raison pour laquelle nous 


la notons différemment ! 


Tl est à noter que l'on peut dériver la relation de récurrence de la fonction J : 
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à partir d'une des formules de récurrence des fonctions de Whittaker : 
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Relation de récurrence pour la fonction d'onde de Bessel modifiée I 


Etant donné que la relation de récurrence de la fonction d'onde de Bessel J est la suivante : 


)- Р 2(о(р-1(р-3)- 8°х°)/, (а, B.x)-(p-1)B*x*J, le, В,х) 
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On en déduit la relation de récurrence pour la fonction d'onde de Bessel modifiée I : 
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Expression de la fonction d'onde de Bessel Y tilde de deuxième espèce en fonction de Whittaker W 


Nous reprenons la formule de construction de la fonction G à partir de la fonction de Whittaker W : 
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Ce qui donne pour la fonction d'onde de Bessel de deuxième espèce Y tilde : 
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En appliquant la propriété miroir de la fonction de Whittaker W, on voit que l'expression est 
parfaitement réel. Et que l'on peut encore réduire la formule littérale : 
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C'est l'expression la plus compacte que l'on puisse trouver ! 


Relation de récurrence des fonctions d'onde de Bessel Y tilde de deuxième espèce 


Qu'à cela ne tienne, regardons avec la fonction d'onde de Bessel Y tilde, quelles seraient les 
formules de récurrence : 
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A partir d'une des formules de récurrence des fonctions de Whittaker W et de l'expression de la 
fonction d'onde de Bessel, il vient : 
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Les deux relations de récurrence pour les fonctions de Whittaker W, s'écrivent : 
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Ce qui donne le même résultat : 
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Lors du cas limite a=0, cette relation de récurrence devient : 
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C'est exactement la relation de récurrence déduite dans le cas de la fonction de Bessel de première 
espèce et comme les fonctions de Bessel de deuxième espèce suivent les mêmes relations de 
récurrence, il vient : 


Expression de la fonction d'onde de Bessel de première espèce en fonction de Whittaker W 


Dans le même ordre d'idée on peut trouver une autre expression de la fonction d'onde de Bessel de 
première espèce en fonctions de Whittaker W, en partant de la construction de la fonction régulière 
de Coulomb de première espèce : 
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D'ou l'expression suivante pour la fonction d'onde de Bessel de première espèce : 
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Cette expression peut encore être modifiée comme suit : 
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Comme suit : 
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Pour ce qui est des fonctions d'onde de Bessel modifiée partant toujours de la définition des 
fonctions d'onde de Bessel de première espèce : 
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En appliquant la relation avec la fonction d'onde de Bessel modifiée I : 
L (e, B, x)= i? J (a, B,—ix) 


On arrive au trois expressions de la fonction d'onde de Bessel modifiée I : 
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Sur le choix judicieux des constantes pour la définition des fonctions d'onde de Bessel de première 


et deuxième espèce 
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Il est important de remarquer que la limite suivante est égale à 1 lorsque a->0: 
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En effet, cette limite va également être étudiée par la suite : 
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D'une part les fonctions d'onde s'écrivent à l'aide des deux constantes comme suit : 
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Et d'autre part du fait де Іа convergence des deux constantes lorsque а->0, rien n'empêche 
d'utiliser une définition alternative des fonctions d'onde comme suit : 
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Puisque LimJ (a, B.x)=J „(B x) et LimY (a, B,x)=Y,(B x) 


Toutefois un autre comportement asymptotique lorsque a est très grand permet de justifier du bon 
choix initial des constantes. En effet intéressons-nous à l'approximation de la fonction d'onde de 
Bessel dans ce cas, et plus particulièrement le développement au premier ordre pour a très grand, 
lorsque x->0 : 
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On voit que ce développement est indépendant de a et que tout autre choix de constante aurait 
conduit à un développement dépendant de a, soit tendant vers O soit l'infini. 


II en est de méme pour la fonction d'onde de Bessel de deuxième espèce, toujours pour a très 
grand, le développement à l'approche de la singularité x=0, et en supposant que p/2 n'est pas un 
entier, ce qui en l'espèce ne réduit pas vraiment la portée de notre remarque, car la singularité du 
sinus disparaît dans le calcul ! 
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Là encore le développement est indépendant de a. 


Pour ce qui du développement de ces deux fonctions d'onde lorsque a et x très grand, il vient : 
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Et les deux fonctions s'annulent bien lorsque x tend vers l'infini, quelque soit la valeur de œ aussi 
grande soit-elle. 


Formules de connexion des fonctions d'onde de Bessel de première et deuxième espèce 


On connaît les formules de connexion des fonctions traditionnelles de Bessel de première et 
deuxième espèce : 
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On connait également les formules de connexion des fonctions de Whittaker M et W, à savoir : 
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Avec les paramètres idoines, liés aux fonctions d'onde de Bessel, il vient : 
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Appliquons ceci aux expressions des fonctions d'onde Bessel de première et deuxième espèce : 
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La formule de connexion inverse s'écrit : 
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Avant de poursuivre les calculs, on peut s'assurer à l'aide de Mathematica, par exemple que lorsque 
a->0, alors l'expression suivante tend vers 2z. Par la connaissance du comportement asymptotique 
de la fonction Gamma, on peut également calculer cette limite : 
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La formule inverse donne également en limite ceci : 
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Cela permet de s'assurer que l'on retrouve bien la formule de connexion des fonctions de Bessel. 
Les formules de connexion sont donc les suivantes : 
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Équivalence des constantes à la limite a->0 


Comme nous avons l'identité des deux constantes suivantes à la limite a->0 : 
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On peut donc construire une équivalence des constantes pour toute valeur a, constantes qui seront 


égales au passage à la limite 0, comme suit : 
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Ainsi donc des deux définitions des fonctions de Bessel de première espèce : 
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Et de deuxième espèce : 
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On peut tout à fait substituer une autre constante pour les fonctions de deuxième espèce, comme 
Suit : 
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On peut même faire de même avec la fonction d'onde de Bessel de première espèce : 
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Avec la relation entre les deux fonctions de deuxième espèce : 
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Et la relation entre les deux fonctions d'onde de Bessel de première espèce : 
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Les formules de connexion initiales s'écrivent : 
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Avec les secondes fonctions de deuxième espèce, elles deviennent : 
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Ces formules de connexion , on peut encore les écrire comme suit : 
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Par ailleurs on peut montrer l'égalité suivante : 
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Cela permet de normaliser les formules de connexion comme suit pour les premières fonctions de 
deuxième espèce : 
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зра) р) е d Ni GIE "Au T ENT 


B [ 52 NP 
Ze * Ë) ES E z 
Sin(pz )/ Jo. B, x)= = - | - SS LE d (co Capa) елэ) 
d SE H Р | d 

a 


Et pour les secondes fonctions de deuxiéme espéce : 


a 


„(ву 
" p 2ле «| 
Sin(pr F (a,B,x)=| Соз(рл)+е 2 J (a, B,x) 2 J le. Bal 
: | ПР, || pH! D | 
2 4a 2 4a 


з}. à "E = 
Sarre ора} Ï (а, Bach x Y 


Mais entre les deux secondes fonctions d'onde de premiére et deuxiéme espéce, la formule de 
connexion devient : 


sc E) 2 
f à Да H x P np < 
Sin(pz )Y (a, В,х) = Cos(pr)+e * Jo. B,x)-2e COR PIE 
r ptl p r p+l Ss 2a 
2 4a 2 4a 


{жы г £) , 
i E z 
Susp Rent EBET ИИ 


Construction des fonctions d'onde de Hankel (Bessel de troisième et quatrième espèce) 


Résumons la situation en ayant développé les expressions des fonctions d'onde de Bessel de premier 
et deuxième espèce J et Y « tildé » ou « double tildé » . En terme de limite nous savons que : 


Lim J (a,B,x)=J,(Bx) et LimŸ,(a,B,x)=Y,(Bx) et LimY,(a,B.x)=Y,(B x) 


a—0 
Et d'autre part les fonctions usuelles de Hankel H(1) et H(2) sont définies notamment par la formule 
de connexion : H,(B x)= J, (B x)+iY,(B x) et HUP x)= J,(B x)- iY, (B x) 


_ B „pH Be fl 
Де Br] 4 2 А 4e Br] 4 2 
J (a, В,х)= e С Ве n 7 паат "aen = 2 Ве n 7 таай Ba) 
(aj: * [r ptl jp "2 (4a): * |rp p*l p KE 
2 4a 2 4a 
„Ë in?! 2 Ë QD. 
ES 4a RP 4 4a RP 4 2 
Y (a, ,x)= E ^ d - 2 тает "ue SE P U ° 7 Wicker ia 
(qay: * |pPl P" a 42): * |g PP. d 
2 4a 2 4a 


On écrira de méme : pri) _ e (2) » E 
H, (0, В,х)= 7 (о, В,х)+1 (о, B, x) et H; (a, B,x)=J„(a, B,x)-iY (a, В, x) 
Aussi je propose la définition suivante des deux fonctions d'onde de Hankel pour un argument x, 


réel positif, œ, 6 réels de même : 


| P Pt 
nb ( R2\2 4 2 
Hie B,x)= 2e sa| P l B WhittakerW P MER. 
S 4a) Jax [p+ D 4a” 2 
I] —+i — 
2 4a 
p? 2 P s DH 2 
m 2 | 4 
Ha, B,x)= 2e sa| B E WhittakerW P P >ë 
i 4a Мах у. pel 7 4a 2 
2 4a 


— Re(HankelH1FromWhittakerW(a, B, p, x)) 
— -Re(HankelH2FromWhittakerW(a, B, p, x)) 
— Im(HankelH1FromWhittakerW(a, B, p, x)) 
—— Im(HankelH2FromWhittakerW(a, B, p, x)) 


Construction de la fonction d'onde de Bessel modifiée K 


Pour la fonction K on va partir de la définition précédente des fonctions de Hankel : 


2 E 22H 
BCEE 4 : 
не Вх) es E L 71 — Whittaker iË P ix? 
4a) Jax Í p+1 В 4a 2 
EISE 
2 40 
p 2 р ia Pl 2 
-rb 2 4 
HO(o,p,x)-2e ^| Ë | .1 É WhittakerW| -iÊ_,P.,ia x? 
40 Мах ү. PHU = P" 4a 2 
2 40 


Avec les formules de connexion suivante pour les fonctions d'onde : 
T Р” (1 
K (a. B.x)- 7 Re ie ? Н("(о, B,i x) х>0 
л ux 
K (a. B.x)- 7 Re ie ? H; (a, В,-іх) x< 0 


Ces dernières sont tirés des formules de connexion des fonctions de Bessel K dont la fonction 
d'onde est censé tendre lorsque a->0 : 


гт 
K (B х)= ie * HVG Bx) x>0 
л o go 
K (B x)-—7 ie H; (=i B x) x< 0 


Comme K est réelle, cela signifie que la partie imaginaire des deux expressions : 


t НОИ 7 ез REI ix) 
E p (a, B,ix et ER p (a, B,-ix 


est censé tendre vers 0. C'est la raison pour laquelle on prend la partie réelle. 


On arrive donc aux deux définitions équivalentes des fonctions d'onde de Bessel modifiée K : 


P Nd B2\2 1 па 2 
K (а, В,х)= 22е а) B Ве 2 — WhittakerW Un 
2 4a) Jax p+1 B 4a 2 
F e qr 
2 4a 
Fes ent 2 
-T — 2 4 
> К (a, B,x)=Z2e sa| В : Im á Whittaker W РВА 
i 2 do ) Jax 4a 2 


2 
2 4a 


Formules de connexion des fonctions d'onde de Hankel à partir de celles de Bessel de première et 


deuxième espèce Y « double tildées » 


Avec les formules de connexion de la fonction d'onde de Bessel Y « simple tildée » 


zB ap? 
Salsa 8,x)= E (zur {г Dee REL 
8° 8° 


Eis ap? * ad 
Sin(pz y (a, B,x)= „(2 J Спірт) bis 
| SS fae fen 25 


3 лв 


EES ОЛ ш te À | Y (a, B,x) 


2 4a 2 4a 


On en déduit celles avec les fonctions Y « double tildées » : 


[ын г”! p 
Pe.) 2 254 2 M ga 
ic) 
2ле es 
4a 
Sous une forme « normalisée « : 
„(рү 
А 2 2ле «E 
Sin(pz F (a. B, x) они) JZ ,x) 2 J uas Bax) 
dër, jB |{Р+1_,В 
2 4a 2 4a 
„рур 
së 2ле i ~ 
Spr) espa d e. al, еа) 
Ë | ‚В | Ë ‚В | 
Г SEA i 
2 4a 2 4a 
En passant aux fonctions d'onde de Hankel : 
, — nep ези 
H; (a, B,x)=J (о, B. x) iY (z, B, x) et [ (а, B,x)=J Je, f.x)-i Ë (a, В,х)© ` н'а, B,x)-H (a, Bx) 
P lepa) SE 


on obtient : 


elita] ll [н 0) (о, В, х)+ Н? (а, В, x) ; 
[zz | .B 


Sin(pr Juta, Ad Ha, A (os eru "sac Ha.) | | 
| | CC 


En simplifiant les expressions, nous obtenons deux formules assez simples finalement : 


Le passage à la limite redonne bien les formules de connexion connues pour les fonctions de 


Hankel usuelles : 


py 
EI „B l: 2. e 2. 8 s 1 
Lim| — | e *T Hi Г Ї =2л S Lim = 
a | 40 2 4a)\ 2 4a Ee pH f. dp) В 2л 
2 4а 2 4a 


Formules de connexion des fonctions d'onde de Bessel modifiées I et K 


Avec les formules de connexion définissant les fonctions d'onde de Hankel et de Bessel : 


af рї Kz: 
I (a, B,x)=e J (a,B,ix)=e?J (a, В.з) (Ha, х) He, ix) 
1.8 1 e 
dë E d Bak dé | e | B 
=1 0, п) = AH Oa, ix) На, Ва) 7 SSC e" ie "e HY(a, В,іх)+ p e 24 H (a, B,ix) 
2л e ta p 
4a 
2 2 
ru (ena E r e , 
I ,(a,B,x)= E E @ "ey i e "re "e ix) 


кә 


š Ne eg (` af 
qr dk LÀ l+ee "e Ha, Вх). — l+ee "e Ha, В,іх)+ 


] LE ИД 
is теде идеала) 


| 


Continuons ce calcul : 
2 2 2 2 
qu JEME y А p qua JEM 
a a ix Log ; ‚фр l a a 
21 ,(a,B,x)= =: xY I (a, B,x)+e” e х1 (a, B,x)+sin(z plie ШҮ d = P гү 
2ле |Z 2ле | — 
4a 4a 
En injectant les expressions des fonctions de Bessel I et Hankel à l'aide des fonctions de Whittaker 
W, il vient : 
DN „АЁ Pt ; NA p 
4a 2 4 4 
I (a, B,x) B x WhittakerW E d + e WhittakerW ч А 
d orl 4a p 1 P 4a 2 p 1 P 4a 2 
I P Dt: 
2 2 4a 2 2 4a 
p D 23 +1 
(Б? \ | RER 2 
Ha, ,ix)=2e «| P á WhittakerW E Pd 
4a va (ix) p+l . 4a 2 
T —+i— 
2 4a 
= PH т?! 
їрл 4 2 рл 4 2 
É : itae Ba |, 6 7 тает A 
, LM (25-2) : (59 A ` 
4a 2 O O 
L (a,8,x)=d Z, p (ов). E 
4a ax \ 4a T ( ) p; re d 
" sin = - fitr iE Piu 
qeu 
y 4a 
2 M 2 P А = PH 2 j z 2 
4a 2 pz 4 -ipt 4 
= E LB TE [з SE 7 mer -iE Bias) é | = EEN 
а ах \ 4a 2.18 a p pel B" а 
2 2 4a 2 4a 
nef) ma 
4 2 4 
Sin(pz Е A set! "Äech 2 j 
ar sud (ar I de TA 
2 4a 2 O O 
= Т 1 (a, ach — E 
4a Мах 4a 
ma а? 
4 2 4 
+i Соѕ(рл £ ; miei id Lau} m | niet - Z J 
pu ps S 
2 4a 2 4a 


| 


On arrive à une formule également assez simple, à savoir : 


| 


2 2 

D. VE А 

1 (a,8.x)= SE S 
LE p E 
2ле <É) 


p 


-gt— 


1+Cos(pr)e 2 


| 


is Bs) = Sin pr)K (a, pa) 


ech 


2 
WhittakerW ү Paso + 
4a 2 


2 
. 2 
P ах 


| 


Il est clair que cette formule redonne la formule de connexion bien connu des fonctions de Bessel I, 


à savoir : ï (ра) = 1 (93)« supr) (8з) = 1,93) L2) Sin(pa)K 93) KB)" Aral 


Il reste maintenant à donner la formule de connexion inverse concernant les fonctions d'onde de 
Bessel K. Manipulons l'expression de la connexion entre les fonctions d'onde de Bessel I : 


(EN estt d ч 
dk | — 2 ж. Ï (a, B,x)= É SH а i 1+Cos(pre * L ,(0,8,3)-2 Sin(pa)K e f.) 
da <Ë] kD P T 
2ле * pete 
4a 4a 


SCH 

ле *| — 

Or (a H d £) | 2° , 4a 
p: 


NES EECH 
bae | SES à [ . do gal Er fat 


x (P озон 37) иеа [ "eege dpt a 


I ,(a,B,x)=C ДЕ 1] Ü Jed 3e [P Бира) 
T a 
B z T B BA . 
>C, — |1+е 2 Cos(pr) K(a,B,x)-K.,(a,B,x)=—C, — е * Sin(pr)l, (a, fx) 
4a 2 '|4a 


K (a,B.x)= C E E Job de sho ena) 


Cela donne donc finalement une formule également assez simple : 


de | Ge 
K ,(0,6,2)= 2 4a 2 4a 


Nd 
2ле 15 
da 


Cela redonne également la formule de connexion connue de la fonction de Bessel K par passage à 
la limite a->0 : КВ КАВ), 


ech Job de Sc 1, (a,B,x) 


L'approche d'Erdelyi et Bateman dans « HIGHER_TRANSCENDENTAL_FUNCTIONS VOLUME 
IL Chapitre VIII, Functions of the Parabolic Cylinder and of the Paraboloid of Revolution» 


Dans cet ouvrage partant des équations séparées de l'équation de Helmholtz en coordonnées, on 
arrive à l'équation donnée au début (modifiée ou non) : 


vide уназ ag уб)-0 


Il est facile de < normaliser > cette équation par une changement de variable qui va nous indiquer 
que le paramètre unique est de la forme + p. 
4a 


2 22 2:72. 
vereor Gr ye ©з E ке) 


Y x 
2772 2: 2 
e Oye Á af ni 
2 y y 
2 2 2 
Posons =; =l к EES w 
a 


C'est l'équation (8.7.1) qui est donnée dans l'ouvrage (avec hélas une petite erreur de facteur 4 
pour le terme quadratique. Cette dernière équation est assimilable à l'équation de Whittaker 
moyennant le même type de changement de variable réalisé précédemment : 


Gy) [s hei 


du 


En posant y(X)= Z(w) = AU) w=7 © AGDE ua t “pop 0 


Posons u=" K =iT v=iw V"(v)+ É d V(v)=0 
y 


Whinakerm( is, P is 


. V 
Solution |, 
Whitaker i, x r3 


Les formules (8.7.4) et (8.7.5) qui sont ensuite données permettent de construire des solutions à 
valeurs réelles, elles résultent de propriétés « connues » des fonctions de Whittaker et des séries 
hypergéométriques confluentes, à savoir : 


TM T их) E d'W. it, — —i J 


-i- (pH) 
zt + 4 4 
d КСЕ: ud v Duch P d in, ix) 
ptl "SN 
(251и) d d 
2 


QA 
V = 


ee nU, ir L ix =2T(p+1)e” Re W ir, P. ix 
s (E 2 


En d'autres termes la fonction de première espèce se construit à l'aide de la partie réelle : 


T 
i (p+1) 
e^ 


T ge D lau 
* (Azur) 


Et une fonction de deuxiéme espéce linéairement indépendante avec la partie imaginaire : 


T 
i" (p+) 
Ceci) 


Ge — mir bi) 
x dein 


Ces deux résultats sont parfaitement cohérents avec ceux obtenus avec le premier jeu de 
paramètres pour la fonction d'onde de Bessel. 


Expressions des fonctions d'onde de Bessel de première et deuxième espèce dans les notations 
d'Erdelyi 


Pour les fonctions de première espèce : 


(г) л(1+р) р 
J (a, B,x)= J (t, %x)= E rd e” $ reri Ze) 


Г(1+ р) 
B? 2e" —т)? 40 
т=-— X=Va х J (a, ,x)=J„(r,%)= ` Re rater BA) 
4a x e e ) 2 
r E ES 
2 “(= y: (=) 
J (a, sexe é = Im —< WinakerW ic, P 
X {> | 2 
L3 Tit 


2 аР 
= лр р 1 4 2 
Y (о, В, x)= ZUM = 2 rater -if 2. ia 9 
гес. x ` T p*l p a 
2 4a 
= ze PHL 
= z^ Р | 4 2 
Y (a, B, x)= “pe, m P EEN 
гес: x p2tl 5. a 
2 4a 
Il vient : 
P "ym! 
= zy EE 2e"(—-r) 4 EN 
Y (o, В,х)= Y, (r. X) =- < C 1)? Re z WhitkerW[ cie Zu 
Á (251и) 
2 
р іа РЇЇ 
= R n 2e"(—r) 4 s 
Y (a, B,x)=Y,(r, z) =-* C M = Whitaker | «is Pi) 
+ do + =) 


Expressions des fonctions d'onde de Hankel (troisième et quatrième espèce) dans les notations 


! 


d'Erdelyi 


A partir des expressions des fonctions d'onde de Hankel : 


2 E rt 
-nÊ ( B? \2 4 É 
HU(,B,x)- 2e "|| . á WhittakerW| i Ê= P. ig x? 
4a) Jax [p+ P? 4a 2 
Г = +i 
2 40 
p £ pH 
_ B 2: \9 1 4 2 
Ha, „x)=2e 4 B e WhittakerW -iP Lia x? 
4a Jax, Pei pr da 2 
2 40 
Оп obtient : 
el 


Whitaker - іт, — „—i a 


шшен fur 


2 
* (Azur) 
2 


Expressions des fonctions d'onde de Bessel modifiées dans les notations d'Erdelyi 


p 
235 PH 2 
I (a. ,x)= | В e * WhittakerM HS 
r( 4a 2 


1+ p}Vax\4a a 
gef 2 D =j, RH 2 
4a 2 4 
I (a, N В Ве £ WhittakerW Ed aa 
Мах 44a p 1 P 4a 2 
D----i— 
2 2 4a 


+л—— P ir É > 
4a 2 \2 4 2 
I (a, r. Ë) Im < reg! 10 vias] 
a 


2 
p2,.1.8. 
2 2 4a 


Il vient : 


S 1 р iP p 55 
1 (a, B,x)- 1, (c. X)e (tne ê Mers is ) 


T(I р) 
шы 
Lf Bao Ls 3)- E CE RÀ =^ ° ниети, 2, im 
a, B,x)- I, (c, X)» —— (=r)? Ке = = >>> WhittakerW| ir, — ,-i x 
é ý X p l. 2 
D —+—+їт 
25 2 
rl 
S 2 -лт р 4 S 
I (a, B,x)- I (r,x)- E (- r): Im ———— >>>  WhittakerW ir Lx 
x gd X p l. 2 
Г| =+=+іт 
2 2 
Et pour les fonctions d'onde de Bessel K : 
Broin? Id 2 
-л—— 2 4 
К (a, ,x)= ze sa| P : Re is 5 WhittakerW D P а 
4a) Jax p+1 .B 4a’ 2 
pu 
2 4a 
p? š  (p+1)z T 
TI — 2 4 
€ K (a, ,„x)= me 4 В 1 Im £ = WhittakerW pr recede 
4a) Jax p+1 .p 40 2 
r E Supr 
2 4a 
Il vient : 
p (p-1)x 
TET m 4 
K (a, В, х)= K (r, x) = Z* ( г): Ве < Whittaker W ir, P. ix? 
А : X (en ga 2 
D ——-ic 
2 
y ;(p+)z 
TE Zei A 
К (а, B, x)= K (c, x) = Z* ( т): Im — WhittakerW| —ic, P. ix? 
p p X EEN 2 
2 


Le processus limite qui permet de retrouver les différentes fonctions de Bessel devient ici dans les 
notations d'Erdelyi : 


Lim J (a, B.x)- J „(B x) Lim J (s, P \- J (Bx) 


Ne 
tin pes) 08) Lim s. Zus 
fin pa)- 1,3) Lim F [s =x - Y (8x) 
Lim Hj (e, В,х)= HY(B x) иин}: ams 
Lim H (a, B.x)- нова) Lim HE E x) неиз) 
Lim 1 (a. B) 18) Lim (r - Ax|=1,(83) 


Lim К, (a, 8,х)= K „(B x) Lim K ут 


PTE NN 
€ 
T = 
9 
> 
Sarre 
Il 
Be 
= 
S 


Wronskien des fonctions d'onde de Bessel, calculé à partir des fonctions de Coulomb 


Partant de résultats connus des fonctions de Coulomb, il vient : 


WiF n x). GG x) - -1 


N 
= 
\ 
< 
= 
> 
< 
= 
> 
= 


Comme Y (a, P, x) = Y (a, B, х)+ fla, DU (a. B, x) > WU Ј (a, B, x), Y Jo, B, x)}= <Â 


Et l'on retrouve le même résultat que pour le Wronskien traditionnel des fonctions de Bessel de 
première espèce J et de deuxième espèce Y, à savoir : 


WU. (a, B, x) Y ja. В.х))- — 
WU a. B.x) Y, (a. B.x))- < 


Avec la fonction Y « double tildé », il vient : 


ehe M s] ә) 
b ;P Hz L) 4a 


2 4a 2 4a 


Et entre les deux fonctions J et Y « double tildé », on retrouve le premier Wronskien : 


„b - 
= (ү (pH .B')(ps1 P 1 2 
J (a, B.x)=< (Ë) Е JL Ey, Ј Да, В,х) А (a, B, x), Y (а, B, x) = 


Il n'y aurait pas d'intérêt à définir une fonction de Hankel < double tildée > à partir de la fonction 
seconde fonction d'onde de Bessel J. S'il en était ainsi alors le Wronskien entre les deux fonctions 
d'onde de Hankel serait de la même valeur que le Wronskien entre les deux fonctions de Hankel. 
Mais nous aurions de la difficulté à définir facilement les fonctions d'onde de Bessel K. Aussi 
contentons nous de donner le Wronskien entre les deux fonctions d'onde de Hankel telles qu'elles 
ont été définies : 


HY(a,B,x)=J,(a,B,x)+iY (о, В,х) Н0(а,В,х)= 2, (a, В,х)-1 (as x) 


TT 
2 40 2 4a 


Aussi de la préservation de la fonction d'onde de Hankel telle que définie initialement on en tire le 
Wronskien des fonctions d'onde de Bessel modifiée I et et K à savoir : 


0. РҮ 
На, В,іх)= J „Да, В,іх) +10 (а, Bi) => WB os, B,i x),J (a, B.ix) = = == a SS p 
3 да 4g 
Ka pale Endo ts) Fla d tJ (a ps) 
э {К (a, B.x).1, (n. B, x))= ДЕ e te p.i J), (a, Bi ) = Rg Fie e la B,i x),J a. p.i 3) 
E | 
2 T. () А d E je je p 
RE 0) S (n n E ME 
Г Fi Г і 
2 40 2 4a 
E : 
m 2ле "^ p 
ME" mam т)“ 
2 40 2 40 


En définissant une fonction d'onde de Bessel modifiée І < double tildé >, comme suit, on retrouve le 
wronskien connue entre les deux fonctions de Bessel : 


ptt. B pe B 

2 Zai 2 E) p x 
| | la Br) eba ав.) 
4a 


Ï (a, В,х)= d 


d 
x 


Formules de récurrence complémentaires pour les fonctions introduites J, Y et I < double tildées > 
ainsi que la fonction d'onde de Bessel K 


2 2 
А ik HZ E 
Posons : cl + а z 
a 


NENG 
2ле (E 
4а 


Сотте : 


del, dei Ei 
Ï (a. B, x)= — ый ; Е A J (œ, aec E hn ie 
e E " 
4a 
B f oi 
EA | 
= a zi 25 АУ 
В TUNE: ИИ F Je. B,x) 
2 4a 2 4a P| 4a 
(zt {г E 
Ee E Lesen Datei ra 
xe" E) н 
4a 


On peut tirer de la constante introduite les expressions suivantes : 


aue ena D py 
2). 2 40) 02 40) "\4a)_ (2а _ l 


Ae QE a2 V 2 ; SV à ; 
2ле (E dE (p-1) JE DE (p-1) 


La première relation de récurrence sur la fonction d'onde de Bessel J : 


J (a, B,x)= (p -2Y2(p -1Xp -3)- P 
Bx'(p КЗ (р 1) 


M IR 
g^ 
= 
F 
= 
= 
— 

CS 


se transforme comme suit : 


Len АР—2Х(р-1(р-3)-8°%х°)9 ү, Le, (20) ар) 
J (a. B.x) gr (p-3) J Yo. B.) P [ | Jo 3) |Z ‚В,у) 


Се qui est la même formule que pour la fonction d'onde de Bessel У : 


a.p. A. eoo -0- P). fa g „y (р, , (2a ет 
ar Aan = (22) o dä вл) 


La relation de récurrence de la forme < double tildée > est identique à celle de Іа 
fonction d'onde de Bessel J, soit la suivante : 


Bell 


nm. 2Y2(p-3Xp-1) B) lap e 
eat ЖЕ (p y). | 
8° 


Et pour ce qui est de la fonction d'onde de Bessel modifiée I : 


sa Xo-29o-1o-3^ 8), C + a. e) 1 E 
P p P В? р 


Elle devient : 


5  2(p-2Y2(p -1Xp -3) вх?) (p-1)f, (2а 
1,(@,ß,x)= Bp) І, Jo B.x) => É IL J2 Jo. B, x) 


Pour ce qui est de la fonction de Bessel K, étant donné que la fonction est construite 
avec la fonction Whittaker W et selon la même constante en dépendance avec le 
paramètre p, la relation de récurrence est identique à celle de la fonction d'onde de 
Bessel modifiée I, à savoir : 


= s= E) a 


Quelques formules limites peut-être méconnues issues de la définition des fonctions d'onde de 


Bessel 


L'expression de la fonction d'onde de Bessel J, à valeur réelle : 


p 

2 2 fl 2 

J (a, B,x)= | А l e ^*MWhittakerM e P ia x? 
Г(1+ p) Vax do 2” 


entraîne une limite peut-être méconnue, d'après la forme alternative de la traditionnelle fonction de 
Bessel en terme de fonction de Whittaker M : 


in 2PH 
e 


EE 


WhittakerM (0, p,2i B x) 


2 p+1 1 
sgg epe . p+ а 


А Ж _ e 1 1 f А Due g€ p 
itio eniti tle 


2 
WhittakerM sh x? 
ax 4a 2 


p 2 Bp LE 1 ( ) 
= Lim = L hs -i „>, ,iax ——— WhittakerM (0, p,2i B x 
eil 4a) Jax 4a’ 2’ 2 [Bx PE 
B° p 
WhittakerM in € dax L 
< L n - 2р1 x WhittakerM (0 0, p,2i B x) 
a? HOER 


De même les expressions des fonctions d'onde de Hankel entraíne également un autre limite dont je 
ne saurais dire si elle est également connue : 


p ; PH 
„8 в2 ү | Dra 2 
Ha, B,x)= 2e sa| P £ ; WhittakerW Ë ЕЕ. 
4a) Jax [p+] .B 4a 2 
Г| —+i— 
2 4a 


Toujours d'après la forme alternative de la traditionnelle fonction de Hankel : 


E 
Hg x) mg pere x) 
TX 


p p ; 
2 2 2 > 2 r -іл 1 2 
Lim H! (0) (а, В, х)= HË ul (в х)=е Ba — WhittakerW(0, p -2i B x)= Lim p S ° Whittaker W i P ia à 
үл а 


a0 4 2 
S Jax T pr B 
2 4a 

р P Lo E -ia x 


i B 2g 4 Ae ` je 
e Lim - WhittakerW\ 0, p,- -2i x) 
a0 | Ae lax e £) 
Г É +i— 


4 


2 
o Lim £ > Hier i 5 ia ch | (> Whittaker W (0, p,-2i 8 x) 
а T 


aa PH 
a’ T DEL e 
2 4a 


On peut encore légèrement transformer cette limite : 


LE Pa TT GI 
4a І 2 4a 
Lim pu +1 2 z Lim 2np^" 
a? de +i E) 
2 4a 


1 
pl 2 2 Ae p-- 
= Lima 2 T EE WhittakerW p pan. m. B NT WhittakerW (0, p,-2i В x) 
a0 2 4a 4a 2 2 


On peut aussi tirer une limite de l'expression de la fonction d'onde de Bessel modifiée I : 


2\2 PH 2 
1 (a, B,x)= = T >= B e * WhittakerM м Ps 
Г(1+ p}ax\4a 4a° 2 


Toujours à partir de l'expression de la fonction de Bessel I d'après la fonction de Whittaker M : 


_ WhittakerM (0, р,2В х) 


BE эге р)ррх 


1 2 jg P 2 
Lim 1,(a,B,x)=1,(Bx)= WhittakerM (0, p,2 B x) = Lim x E 1 imate | ia) 


2" T(1+ р) |28х 4a 


P 

1 232 j PH 2 ; M 2 

© Lim — a e * WhittakerM РВ x? s (0, p, Bx) 
a0 lax 4a 22Р 28х 


-ї—,‚—,— 


2 
тиде "e P! — Whittak M (0 2 ) 
a ir — 
o Lim p+ =e 4 Ух — c = 
a0 a? (20): 


Il en est de méme avec l'expression de la fonction d'onde de Bessel modifié K, à savoir : 


(р=1)я 
Fg 4 2 
K (o, В,х)= z e sa| B | Ке E = WhittakerW D D sup 
Vax (224 J 4a ” 2 


2 4a 


Ce qui conduit à la limite suivante : 


K, (B х)= (Eater pa x) 


NH 4 2 
Lim K (o, B,x)- К, x)= — WhittakerW (0, p.28 x) = Lim z e <| P l Rel d zx WhittakerW РР раз 
а0 2px Jax [E D | 4a 2 
Г Ыы ты 


WhittakerW (0, p,2 x) 


P 
_ B 2\2 ] ar 2 
© Lim e «E Re S ritieni ia) = 
a 


1 
Vax {жык J2p xz 
4a 


1 
2 P 
o Lin S Re 5 нийет iE 2 ia 3 = 2 ec © WhittakerW(0, p.2 0 x) 
m pt, Л 


4a 2 


1 
p 1 2 (p-r BU? De 
= Lim Rd a 2 E e * WhittakerW| i га y: ia x {à fr xWhittakerW(0, p.28 x) 
a 


Comme les parties imaginaires de ces deux limites sont censé être nulles, on en déduit : 


See Arie | 
E 2 1 4 
Lime * L S WhittakerW PP y = 
dor р+1 Ж 4a 2 
2 4а 


P k M 


p p 2 Ste 
Whittaker W бе т, =— = WhittakerW (0, p,2 x) 
O p+ NT 


B 2 


WhittakerW (0, p,2B x) 


1 
\2B xz 


Z [psl g) et p "i 
Lima ° (ear 4 minen EL Gg ax | É | Je xWhittakerW (0, p,2B x) 
4a 


(04 


Avec le formalisme d'Erdelyi, ces limites deviennent, pour celle issue de la fonction d'onde de 
Bessel J : 


2 B? p d 1 
En аем @ 4 = Ia x -_ = WhittakerM (0, p,2i B х) 
=> ax 


a 2 — 2 bg: 


2 pH 2 WER 
T = p > Lim t? Vio r 2, б) ei = Bx WhittakerM (0, р, 2i B x) 


до, T—+0 


;„ P 


pa 2 cH. 
— Lim т? Mise inti i + E 3 x WhittakerM (0, p,2i x) 

T->+œ0 2p+= 

2 2 


4r 


Et pour celle issue de la fonction d'onde de Bessel modifiée I : 


р 
242 j,PH 2 
yp i B e ` ^ WhittakerM s D dd 2 _ WhittakerM (0, p.2f x) 
4a 2 22 2B x 
oe PH 


2 pH 3 
T -Ê Be — Lim t ? WhittakerM| –іт, E = б) x) -JBx — > WhitiakerM (0, p.28 x) 
4a 2 d p 


T—+0 Ar - 


pu 2 p 
— Lim т? тием ins, iz + © 5 Jx WhittakerM (0, p,2x) 
T—+00 T 2p+= 


2 2 


Pour celle issue des fonctions d'onde de Hankel : 


2 
nÈ -Ë нийет i 3 I | 
S (04 -in— 
Lim IÉ à =e * WhittakerW (0, p,-2i B x) 
a>0 | 4a) Jax e .p | 12D x 


4a 


р Bx 
i WhittakerW Hoan P 
F T -in— 
EH уа — Итт?е" =e * BX Whitaker] (0, p.-2i B x) 
4a "ur a>0 [tus] 8л 


2 
„ WhittakerW| ic, Р EI Ps 
р 2 4 ir 


pu np 
> Lim t ie" eist WhittakerW (0, p,-2i x) 
a0 Е +1 I ) 8л 
Г SCH LIT 


Pour la limite issue de la fonction d'onde K : 
2 Ak 2 
EE e * : BP P, 3 l : 
Lime ^ WhittakerW| i —,—.,ia x^ |= WhittakerW (0, p,2 B x 
e r 4a 2 DS (0, p.28 x) 


2 4a 


2 pu 4 — 2 

T m {БЕ > — Lime ^c? ——<————Ийайег ir Ë 8з) = yp» WhittakerW (0, р,28 x) 
4a AP T4900 [ter 2 4r 227 

2 


pu 2 

— Line^c? ea iae = = X шшде, px) 
уш il p+l sr) 2? 4r) „т 
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Dérivées premières des fonctions d'onde de Bessel de premier et deuxième espèce 


Dérivée première de la fonction d'onde de Bessel de première espèce J 


Toujours dans N.Y Vilenkhin < Special Functions and the Theory of Group Representations », 
« Chapitre 8 Whittaker Functions », page 416 et suivantes, il est donné les formules de dérivées 
première de la fonction de Whittaker M : 


ôWhittakerM(v, H, z) = (2u + 1} 


= WhittakerM (v ,u,z)+ ( 


2u+1+2v)\2u+1-2v) 


WhittakerM (v, u+1, z) 


2 i2 [(3+ p) 


д; 22(2и+1) 8(2u+1Y (и +1) 
2 
ôWhittakerM(v, u,z) _ Qu -1) M WhittakerM (v, u, z)+ 2 uWhittakerM (v, u —1, z) 
z 2zQ2u - 1) 
De pl TEF E last) Б j 
e plus rena Е E . Donc : 
Hp у ke 
a ? i? T(1+p) 
2 
ШОШ 2) 
=p P 2 
a = 6А 2ia 5 sie Е Еа 
X Z 23 “rü p) y x Hei 
y ia x? 
B? 
Whi M =Â 
ai le, Bx) 2p | OWhittaker. (Z. >] (a, B, x) B? 
a S| a ra 2ia P. 2 yoiax PE 
ü 2 j ? I(1+ р) 2 
y-iax? 
x 2 
Or OWhitiakerM(v. u, y) _ Qu +1) —2yv WhittakerM (v, и, у)+ Qu LEV Ouen) ana S) 
Oy 2y(2u +1) 8(2и +1) (u +1) 
2 2(p +1)" - gx ; B Pay 
" omiitaer ZZ а тети 
u=? 2yv P => SS = Bt 
2 2 ду (p+1) + B +2 
|н 40^ WhittakerM PTS йал? 
y=iax 4(p+1) (p +2) dia” 2 
2.2 2 
2(p+1) =B x наге SE P dax |+ 
Ме 4i a х 2(р+1) dur 
Ес `> 2 2i a 4 
aJ (a, B,x) _ pp нр (pada B 
e Óx 3n s Tix p) t д0? WhittakerM p? Hope 
4(p +1) (p+ 2) 4ia 2 
_ J (a, B.x) 
x 
2(p+1f -pe J (abs) J (a. fen) 
( ) 2(p+1) x x 
0J (a, B,x 
SEE ôx = S â 3549 „ MWhittakerM BU АЕ. 
„Макы Sp Aia" 2 
2B 2(р+1) 3+p 3+p x 


Soit le résultat remarquable : 


ôJ Да, В,х) _ 2p(p+1)- Bx? 4a*(p+1) + B* 
àv pe DIVI дои ei 


Lorsque a=0, on retrouve bien la formule de la dérivée première des fonctions de Bessel : 


0J,(0,8.x) _ 2p(p+1)- ^x mE 
Be ^ abe Pru pro) 
Fonction de Bessel J 
aJ J 
LD LE gp (ва) Ais (Ba), ala) 
x x B x 
__ Вх a (Bx) J (Bx) f 
= J, (B x De y tx) ex =P x pa 691,202) 
êJ (B x) _ p - Bx m Bx 
=> ES -2 2(р +1) J (Bx) Ap +1) 09) 
êJ (Bx) _ 2p(p+1)- Bx? В°х 
ЕТ 


Avec Іа première formule de dérivée il vient : 


Or ôWhittakerM(v. m.y) _ Qu- 1) eii WhittakerM (v, H, y)+ 2 uWhittakerM (v, n Ly) 


ôy 2y(2u-1) 
2 21.2 
yoiax v> É u=7 уу Et 
2 
тийем A ») 2 
=> = = 2(p- у SEA WhittakerM PD ш + pWhittakerM £ E" ‚Тох? 
ду Aia x’ (p-1) 4ia 2 4ia' 2 
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Z 2 
2(р-1) =P WhittakerM В ,P тах? + 
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J (a, ,x) 
x 
2(p-1} - Bx? J (a. Вх) Jo fo) 
( ) 2(p-1) x x 
0J,(a, B.x) _ m 
dnm ES = А Е WhittakerM ° o а? 
xB 2-2) gr dia 2 
2 pipi 
2(p 1) vja r(p-1) x 
oJ ‚8, 2.0 5d _ 2 
Ps Bx) Вх? -2р(р-1) (зру АВ р (ов) 
Ôx 2x(p-1) ^" 2(p -1) 


II s'avère que c'est exactement la formule de la dérivée pour la fonction de Bessel J usuelle : 


p BJ, (Вх) 225 H AECH, (ва)- л, (в) 
=, (69 3 2 SU, 9), (8 x))=> АӨ) Рз), Aal ZL 
Е ECH, (в), LE EE c.q. f d 


Dérivée première de la fonction d'onde de Bessel deuxième espèce 


Toujours dans N.Y.Vilenkhin < Special Functions and the Theory of Group Representations », 
Chapitre 8 Whittaker Functions, page 416 et suivantes, il est donné les formules de dérivées 
première de la fonction de Whittaker M : 


. 2 
ôWhittakerW (v, u, y) Е Qu _ 1) = WhittakerW (v, u,z)- eus WhittakerW (v, u -1, z) 
ôy 2y(2u -1) 2(2u -1) 
OWhittakerW(v,u,y) _ (2u+1Y -2yv WhitiakerW(v, u,2)- (2u+1-2v) 


Whittaker W (v, u +1, z) 


ôy = 2y(2u+1) 2(2и+1) 


Soit dans la fonction d'onde de Bessel de deuxième espèce : 


p- SIDD VE D E prem £, P triax `) 
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л p^ x 
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у=іа х 
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ЖОШ E 


pA | | 
ato) — Re Ed EE) ра 
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Continuons les calculs : 
Ӯ ( B ) p z ome Р Ру A B ) 
p% P >X) _ „ (4a)? Re (Т р+1 B dia ia Va, B,x 
Ox л p^ 2 4a Oy x 
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"b; 5 i 
à i 2( тр 4 ы = E 2 E 
` КЕ (p-1) ЕВ L 2 (4a) z Re (+i) T РЕ =P WhinakerW |v, Z 2 гах) 
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Ox 
Lorsque a est nul, on retrouve bien la formule usuelle pour la dérivée premiére de la fonction de 


Bessel Y : 
2 A BY, (B x)— p P Y, (8 x)= 0-10), )- Y, (8х) 
m (v (Bx) Y, „(В x)= Z -p EEO Dale) 


Bx) Вх? -2р(р-1) p? 
= ais Y (B x) de r2») c.q.f d 


Avec la deuxième formule sur la dérivée première de la fonction de Whittaker, il vient : 


yoiax e u =P à SE 
4ia 2 2 


5 2 
ôWhittakerW (v, u, у) = (2u+1) -2»v WhittakerW (v, u,z)- Qu V 20) ушел, u +1,z) 
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=> 


(a, B x) 2 p+l е 
=> = Вх (40) 2 p+ p*3, B | : | 
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Ў (а.х) 
(a. B.x) 
x 
ôY,(a.B.x) _ 2p(p+1)- £^ > Bx ӯ 
=> ER = 2x(p+ 1) TLA e ee ,x) 


Il s'avère que c'est exactement la formule de la dérivée pour la fonction de Bessel Y usuelle : 


Fonction de Bessel Y 


OY, T 
= АВА) By (ва) y (p) ШЕШУ (84). y, (x) 
x X px 


= аз) PE EOD (p DEE P Gr (азуу, (px) 


дү (B x) E >x 
KE LE ne ra 


dY (Bx) _ 2p(p+1)- ^x p^ 
=> a 2x(p +1) Y (8x)- 2 SR c.q. f d 


Autres expressions des dérivées premières et secondes des fonctions d'onde de Bessel J , Y tilde et 
de Bessel modifiées I et K 


Ces calculs ont été réalisées à l'aide de Mathematica et redonnent les mémes expressions que 
précédemment, et je compléte avec les expressions des dérivées secondes. A partir de l'expression 
de la fonction d'onde de Bessel J en fonction de Whittaker M : 


e 

4 JP 2 

J (a, a) E Е) 
а 


On obtient les dérivées premières et secondes, comme suit : 


235 2 
B (-2a + 2i xa? «if? WhittakerM 292 P ye + 
dJ (a, B,x) 4a -ir +2 4a 2 
Ë Е Кее * 
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A partir d'une des expressions de la fonction d'onde de Bessel Y tilde à l'aide des fonctions de 
Whittaker W : 


P 

~ 2\ 2 i P 2 2 

Y (a, B, x)=- І P Rese *Г DEL ы WhittakerW =; PS B r 
лхуа \ 4a 2 4a 4a 2 


On obtient les dérivées premières et secondes, comme suit : 
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À partir de l'expression de la fonction d'onde de Bessel I en fonction de Whittaker M 


SECH В? p 
=> Z — Ree * WhittakerM| -i 
Г(1+ pix la 


i——,—,-ia x? 
4a 2 
On obtient les dérivées premières et secondes, comme suit 
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4a 2 


W: 
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On obtient les dérivées premières et secondes, comme suit 


2 
p Их (2120? -2a -iB* WhittakerW EP LE = 
dK а, В,х) ле“ (Z) r e 4 4a” 2 
e 
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r E - 4a WhittakerW [a S, Pia x? 
2 £) 4 
2 
2(R2 , 5% | ‚В p+2 . 3 
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À partir d'une des expressions de la fonction d'onde de Bessel K à l'aide des fonctions de Whittaker 


FE Xx 
4a 2 


2 
L iar) 


Simplification des dérivées premières et secondes des fonctions d'onde de Bessel J, Y tilde, I et K 
Par identification de l'expression de la fonction d'onde de Bessel J : 


P: 
2 


Z) 4 

4 _ LP 2 

J (0.8. x)= 2 Bei e ^WhittakerM hope 
Г(1+ ph Va 4a 2 


On obtient pour la dérivée seconde : 


_ (1+2р) f? +2(1+ p p(1-p)+x*a?) 
2x°(1+ р) 


d'I os fer) Ae pf a" egt 
d? 21+ p)P J „ala. B.) 


J (a. B.x) 


En jouant sur le respect de l'équation différentielle des fonctions d'onde de Bessel, on arrive à 
l'expression suivante de la dérivée première : 


d J (a. В,х) 1 CE 
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dx? x dx 
E" 1 dJ (a. Bal d J, (œ. B.x) ax р? хв? e 
x dx dx x | | 
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Par identification de l'expression de la fonction d'onde de Bessel modifié I : 
4a ER B p 
I (a, В,х)= === Re e * WhittakerM|\ -i —,— ia x? 
À f ) T(1 p)va | | 4a 2 ) 


On obtient pour la dérivée seconde : 


(1+2p)x*8"-2(1+ ppl- p)+x*a?) 
2x?(1+ p) 


d1 lepa) _ A+ pa? + B' 
dx? ж 2(1+ p)£ 1 Joch 


1,(a. Bal 


En jouant sur le respect de l'équation différentielle des fonctions d'onde de Bessel, on arrive d 
l'expression suivante de la dérivée première : 


d'1,(a.B,x) 1 dI (o, B, x) 
dx? x dx 

Ex dl (z, B, x) = 2p p)- x^g? 
dx 2x(1+ p) 


(es -2i p |1,(a.p.s)=0 
x 


4(1+ pf a? + 8* 
2(1+ р)8° 
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Par identification de l'expression de la fonction d'onde de Bessel Y tilde : 


p 
= 212 Py! 2 2 
Y (a, B,x)=- NE ES Rei e IT peut yip WhittakerW E Saa 
л xa V 4a 2 4a 4a” 2 


On obtient pour la dérivée seconde l'expression suivante : 


d Y (o, Bx) _ В? E (1+2p)x28?+2(1+ ppl- р)+ a°] ç 
dé 21+ DE (een) 2(1+ р)х? Y (a. Bal 


En jouant sur le respect de l'équation différentielle des fonctions d'onde de Bessel, on arrive à 
l'expression suivante de la dérivée premiere : 


d Te. Bux) 1 deld ans - P Pg k (a, B, x)=0 
dx x dx 
Sade) Pep) atop FU f 
x dx dx? x? 
dY (a, Bal B à 2р(1+ p)- x? f? 
EE DEUS des pex 2(1+ р)х TEES 


Par identification de l'expression de la fonction d'onde de Bessel modifiée K : 


B? E 
77да 22 m 2 
К (a, DE P Re =< > MWhittakerW ;P P: а 
хуа (4a p+1 .B 4a 2 
D =— +i 
2 4a 


On obtient pour la dérivée seconde l'expression suivante : 


d’K (a, b,x) _ 4a" (1+ р) 4 f* 


(1- 2p)? f? -2(1+ platt p}+ vie" 
dé ` XI p)P Sech ëch Kj. f.) 


201+ p)? 


En jouant sur le respect de l'équation différentielle des fonctions d'onde de Bessel modifiée, on 
arrive à l'expression suivante de la dérivée premiere : 


d'K (a. Bx) 1 => — ech (a, B,x)=0 


dx? x dx 
ак „(a, B, x) _ MN ?(1+ p) + p* x' B? +2p(1+ p) 
Son 2(1+ р)8° K alc; Bax 2(1+ p)x E 


Voici d'autres formules pour les dérivées des fonctions d'onde de Bessel 


dJ (a. B,x) _ B*x* -2p(p-1) 


1 Be 2, еВ) 


Ox R 2x(p -1) 
а°7,(0,В,х) 1 dJ a. B.x) _ a*x* -p +p? , PE 
dx? x dx x? S 
_ 44,0. B.x). Xp-1p(p+1)-a*x*)+(1-2p)B*x? , (o В,х)- В? J 
dx? 2x'(p-1) pU 2(p-1) ? 
dY,(a.B.x) _ В?х? -2p(p -1) Y,(a.B.x) | 4a (p- 1) + B° d b в E 
dx 2(p-1) x Acn" жк 
d’? («,В,х) _ 1 а, (а, В,х) a*x*=p*+ B*x? Y (a, B, x) 
dx? x dx x? AU 


_ 1, Bx) _ 2(p-1Y\plp+1)-a*x')+(1-2p)B*x? ç ( Bax) pP 


dx? 2(p _ 1)x? 


dI (a. Bal — xp +2p(p—1) 


B? 
ES = SEET го, Вх) x> yes B.) 
d'I(c,B,x) 1 dl,(a,B,x) завеса, ss SA 
dx? x dx 
d'I(o.B.x) _ Xp- 1(p( (p+1)- a? ios B? 
> dx? ü SE 1)x? 1,(a. B.) Xp-1) ^? 
dK ,(a.B.x) _ х2? +2р(р-1) „ (aua s) 
dx 2(p —1)x д 2(p-1) EROS 
d'K (a. B.x)  14K,(a,B,x) ax! =p" – В?х? = K (a, B.x) 
dx? x dx x? E SH 
_ PK,(œ.B.x) _ (p-Wplp+1)-a*x*)+(2p-02B" gi Ex 
dx? 2(p-1)x? ИЕ ” 


Pour les fonctions dítes « double tildées » : 


Sp 
А (£u (25 e) S i 
En définissant 2E + ? 20 : 2 20 , et sachant que Teac (f (o, B,x) 
4a ES p p Pl 4g |? 
2ле (E = 
rap o) 
©) 
TL 4a 
Сотте : 
D D 
20 ваар сы). p' 
c. | B В' c( E) Fete y 
E 4a P| 4a 


ó 6 T (Z) 

p Ae _ p’ p+ Aer i B*+4a°(p+1) 
ë 4£ 8° +402(р+1) " a 8° 

PT 4a Р 


Il vient : 


11 


„(a.B.x) _ Вх? -2р(р-1) 5 B^ +40? (р-1) 
Ox = 2x(p-1) 7а, Boxen 


di (a, В,х) e 2p(1+ p)-x?B? = pos 
dx g 2x(1+ p) Jala Bir) =x d 


d'Ï (a, p, x) 2(р-1)(р(р+1)-?х* +(1- ) 


» )+(1-2p)8*x? = f^ «4a? (p-1y = 
dx? = 2х°(р-1) J (a B x)- (p -18? J, Jo p x) 
427 (e, b,x) _ (I+ 2p) 202 +2(1+ pXp(1 р)+ a?) S B? = 
dx? => 2x?(1+ р) J (a. ëch alle Bal 
dY (a, pa Bx? -2p(p -1) Y (a. B,x) B 3 
dx = EE x cies en) 
dY,(a.B,x) _ 2p(1+ р) хв? 5 psa (ps1) 5 
dx 2(+р)х ло 201+ p)g? faa Bx) 
4 Y o. B.x) _ (p-VNp(p+1)-a*x')+(1-2p)B"x" 3 _ B $ 
dx? ü 2(p-1)x? Tres B.) Xp 89) 
d"Y (a, B, x) " ((1+ 2p)}x° 8° +2(1+ plett р)+х“а?))= GER B^ +4a*(p +1) 
dx? 2(1+ p)x? OLD 


di (a, pun) _ xf? «2p(p-1) B'+4a*(p-1f = 


di 2x(p-1) ech: gece elei 
cx un NO а,в.) Е 
1, с _ Xp-)p(p+ =” - E na ESCH, (a. p.) 
j _ (1+2p)x*B? Am ая + doge В.) 


Formules de continuation analytique pour les arguments négatifs 


À l'aide des développements limités autour de zéro, on voit rapidement que l'on peut établir une 
formule de continuation analytique comme suit : 


J (a, B.xe"7)- e" J (a, B, x) 


Il en est de méme avec la fonctions d'onde de Bessel modifiée I : y (a B хе"^ )= eP t] (о, B, x) 
p H H D H H 


On peut établir également ces deux formules à l'aide de la formule de continuation analytique de la 
fonction de Whittaker M (voir NIST Handbook of Mathematical Functions, formules 13.14.10 et 
13.14.11) comme suit : 


WhittakerM (к, ui, ze” )= +ie*" WhittakerM (- , u,z) 
WhittakerM (к, И, ве") (-1)" &"WpittakerM (к, u, z)=e""e*" ""WhittakerM(s,u,z) тем 


Р 
, 2 1 2\2 р 2H 2 TESTI 3 
Avec l'expression уйг В e "T WhittakerM UE Pl , il vient . 
Г(1+ p) 4a) Vax do 2 
2 2 
x хе!"" => WhittakerM m m xig" |= (- 1)" e" WhittakerM PUN x? 
4a 2 4a 2 


inm 
eu 


cela donne ` J (a, B.xe"*)- C1)" — J (a, B,x)=e"*"" J (a, В,х)` 
e 


Il en est de méme avec la formule pour les fonctions d'onde de Bessel modifiée I : 


P p 
1 2 2 j PH 2 1 2 2 oe Bil H 
1 (o, B. x) - ———— —- В e * WhittakerM Р ar == m А e *WhittakerM {Б Р. 
g r( 4a” 2 ER 4 


lux 
2 | 


1+ pax да 4a 


> le. B, хе") е" (a, b,x) 


a 


1+ pax 


a 


P 


La formule de continuation analytique de la fonction de Whittaker W est la suivante (voir NIST 
Handbook of Mathematical Functions, formule 13.14.12) : 


(=1)""27z i Sin(2mur) 


WhittakerW (к, u, ze?!" ) Geer 
É - u- | +2u)Sin(2ux) 


WhittakerM (x, u,z)+(-1)"e ?""WhittakerW(K,u,z) meN 


Voyons si ces formules sont également valides avec la représentation de ces deux fonctions à l'aide 
des fonctions de Whittaker W. Et pour cela appliquons la formule comme suit : 


(1) 2x i Sin(mpr) 


d = d Si (i+ p)Sin[pz) 


2 2 2 
x 2 xe" => WhittakerW De xg" | = WhittakerM PR x «Clem "WhittakerW AP x 
4a 2 4a 2 4a 2 


2 
WhittakerW 2 p dg xdg 
4a 2 , 


2 
dp _;P 
2 4a 


2 
WhittakerW Apnd x? 
4a 2 


2 
Гр. 
2 4a 


(-1)""2xiSin(mpr) 


qub nt 


2 
WhittakerM P AL-Uemg 
4a 2 


2 4a й 
2 2 
De ph Lit VE i, | 7 A s - 
6 : Sin ЕЕЕ 2 Sin SE 
2 4a 2 4a 
2 2 
WhittakerW AP Pu xe | (-1)"2isin(mpz Sin Pg YR WhittakerW Е А P ja x? 
4a 2 2 4a ) Pp., EE 4a 2’ 
- = - WhittakerM| AË ‚ах? |+(-1)" e" 7 
р bep 8 Г(Ї+ p)Sin(pz) 4a 2 r Ip B 
2 4a 2 4a 
= 2 2 
WhittakerW TA pem А leck Sin(mpz )Sin Re | WhittakerW JE Bas x 
da 2 2 a | B p. ety m do 2 
; - ri pni) WhittakerM Love ia x? tem | 7 
Г Ee FPE Г e 
à 4a à 4a 


Pour cela prenons les fonctions d'onde de Bessel J comme suit : 


2 а В? 
(Ë) | e ta Р 2 = EHE 2 
4a 4 4 
J (a, B,x)= £ 5 WhittakerW y Puy + S x WhittakerW GE E ua 
хуа dei P do 2 det P do 2 
2 4а 2 4а 


Sachant дие : 


Ve B p s Dr 
e *MWhittakerM| -i —,—,iax. |=е *WhittakerM| i——,—-iax |= 
4a 2’ 4a 2 


É gna Tr p) 


Ke 
H 4a 


2 2 
Whittaker W Ë „Р ja x? +p) WhittakerW bus 
Ae 2 4a 2 


2 
Г pA C 
2 4a 


Cela donne : 


e 2 
2 a C1)" 2i Sin(mpz )Sin p Dyk. EE 
B 7 LE 2 4a J iL -m “> P. 5 
CS е “а гіх Iw ) e ?e *WhittakerM SE + 
j іпірл j 
J (a, p,xe"" )= E | j 2 +e J (a, В,х) 
xe™ Va T 
LU" 2i Sin(mpz )Sin z wf У 
2 4a РЕ j;2H p p 2 
+ - e ?e *MWhittakerM| i —,—,-ia x 
Г(1+ p)Sin(pz) 4a 2 
2\5 _ BË p 2 
2 -r> -in PE 
изип] Ë | e g + rien -if аз E ў E 2 
J (a, В,хе"")= us M e 2 Sin Е: re +e”? Sin TARTA 
Sin(pa)x/aT(1+ p) 4a 2 4a 
te" (a, B,x) 
sE 1 2 l 2 
ЖАД 4a гі _ sip P й 
J (a, В, хе" )- 2i Sin(mpz je J (a, Che 2 Sin P i YE d +e ° Sin z же Tte (a, B, x) 
? Sin(pz) б 2 4a 4a I 
NA ze PH pet 2 aH — 2 - P -ij m i ў mrb 
Comme e "l LOST "(eta se ше "ais 
a a 


> J (a, B, xe"*)= 2i Sin(mpz J (о, В,х)+е ""*J (a, B,x)= e"" J (a, Bal 


Il en est de méme avec l'expression de la fonction d'onde de Bessel modifiée I : 


2 


4а 


2 
vue! i Bp ia 3 4—————M 
n" | 4a” 2 d B 


2 


Passons maintenant à la fonction d'onde de Bessel Y tilde, telle que défini comme suit : 


2 
Whittaker W шы, x? 
4a 2 


ia x? 
2 ) 


ах) 


р apt 
zi 1 Ee 2 1 2 LE 2 
Y (a, В,х)= P pP Y Egg e WhittakeWw| -iL_ 2. 
z xVa | 4a 2 2 4a)\|2 2 4a [ptt T 4a 
2 4a 
Soit : 
i PH 2 
4 
, Я „2 = Whittaker -i£ 2 ias - 
2 +1 ` a 
E pcd E: (25-2. 
> 2 2 4a 2 2 4c J) 4a 2 4a 
Y (a, B,x)= - " 
2iz хуа SE 2 
E = тикети | +i Ë L 
2+1; z 
2 40 


Il vient : 


mH a. Q: Š 1-р ó p 
2 (ч) ime Doc. "A j 
r| 24 | Di r 24 | D py 1 Ç M q Lg E ment. Pi „idx 
2 2 4а) (2 2 4a \4a Г(1+ p)Sin[ pz) 4a 2 
V imz | _ ; imz 2 
Tee J. E "2i Sin[mpz )Sin Fu N 
P 2 40 Pt Dag 2 
Гї+ 5i I ) e Ze? ZZ Siar 
p)Sin pr a 
+e "Y (o, B,x) 


e P Y (o, B,x)+ 


B rg pu 2 

: NT Ph LN SE 
ac) J T| IE L Ve *MWhitakerM -i E iax |le ° Sin z +I [E 

> Y (a, B,xe )= | É 2 E) [ 2 4a Sin(mpz) 4a 4a? 2° 2 4a 


Jr > sen) ШҮ a ет, 
Ag 2 4a 
+1 2 
lÉ Е EN S A aftu 
(T imz -impr үу 2 H 40 2 2 4a Sin mp T ; ч 
> T (Ве )=e p Y (a. 8,x)+ Дт: SE NI l- 2 
Л p -e "2 Sin =L EB 
m 2 4a 


, pH 2 , pH 2 B „p „p B 

ig 1- ir = f= -7 = -r 

Comme e ? Sin| z He -e ° Sin E -ie * + Cos(prk da | = je а Соз(рл)+е ? 
2 4a 2 4a 


2 H 
PP up - P 
2 2 dajl2 2 4a)Cos(pr)+e * 


SERES 


1 
= Ў (a, p,xe"" E e "Y (o, B,x)+2i Sin(mpr) J (a, B,x) 


On retrouve bien la formule de la continuation analytique de la fonction de Bessel lorsque a tend 


vers 0, à savoir : 
2 2 
[zaid LEE 74 


| =15 A xe"*]- ey (B х)+ 2i Sin(mpr) 


Lin 2 2 4 4a cal a (nd 
dÉ e ta 
4a 


La formule de continuation analytique est donc : 


2 
{ын EE 
F (a, B,xe"*)=e "rï (a,8,x)+2i Sin(mpz ) o | u Cos(pr)+e 
| | QE Sin(pz) 
gi e 4a 


а 


J (a. B,x) 


) 


| 


Voyons maintenant ce que cela donne avec la fonction d'onde de Bessel modifiée K : 


"d 2 Р El 2 _ B 3 р dei, 
4a 2 4 4a 2 
K (a, ,x) =Z p £ zx WhittakerW yes SEZEN Е 

xia \ 4a p+1 D 4a 2 lixa 4a -iz 

T| =——+i— ó 

2 4a x | 
е E 

2 


Nous savons que : 


i D 


4 
= = Whittaker d I 


Eg UE 4 


Bp 


а 


EP 


а 


SNE 
IG X 


} 


. 2 
IG X 


2 i е, 1- 2 2 
пайа |17 а) (у) озне sel Я WhittakerW E L а 
| = ; = x por) s mao À Äech (Kg = - 
п 2 
hp. B° p Siapa d nbus 
2 4a 2 4a 
2 m z: LS 2 2 
meter caen (-1) Seel Zuel з тийн I 2- 
; - S GE Е Ë ia ch (eg ; 
Jan 8 TU: p)Sin(pz) 4 dp e 
2 4a 2 4a 
On a donc : 
2 
Eater Er | R : 
LE je M I | | i P mine Pia + 
| | 4a 2 + pJòm\ pa a 
K Да, В, хе") е" (a, Bh Z< E 5 
2i xe" Ja | 4a Er. ере, B 
(-1) 2i Sin(mpr )Sin EE KEE j 
2 e le *WhittakerM 4T 
Г(1+ p}Sin(px) 2 
B о а? B Иш 1-р p 
p Sil II g^ imi Pia e wll 
: К -л— O (0 
= K (a, fae") e^" K (a, 8,3) "n^ S Ge WI ) pH 2 
in pz tpa sr = 
Р 2 -e ?Sin EE 
2 4 
Bo NS 2 NS 2 
; 3 ^ Sin(mpr) шоны e LP e Fali E 
> K la,B,xe"" |=е "К (a, B,x)-re * I (a,B,xKe ? Sinr—-irx —|-e ? Sin x —-+ir — 
/ B ) / p ) OD À p 2 4a 2 da 
рі 2 pd 2 p p 
ir — ]- -in—— ]- T— =л— 
Comme е ? Sin SC -e "Am pb -ie **31+Cos(pr)e 2 
à 4a 2 4a 
NA 
1+ Соз(рлђе 2 


> «(ЭЛЛЕ e""K (a, В,х)- ix Sin(mpz) 1 (a. B,x) 


Sin[pz) 


La formule de continuation analytique est donc : 


=й 
1+ Соз(рт)е 2 
Sin(pr) 


Et l'on retrouve la formule de continuation analytique de la fonction de Bessel K lorsque a tend 
Sin(mpz) 


vers Ü : K (в хе") е" (B x)-in = 1,(Bx)- 
P ? Sin[pz) ” 


K la, B xe"*)= e" K (a,B,x)-ix Sin(mpr) 1 (a, B,x) 


Pour la fonction d'onde de Bessel Y « double tildée », la formule de continuation analytique est des 
plus simples : 


Ë 
Соз(рл)+е ” 


Sin(pr) Alstad 


F (a, p, xe" )= e" (a, В,х)+ 2i Sin(mpx) 


Pour ce qui est de la continuation analytique des fonctions d'onde de Hankel, étant donné leurs 
formules de construction : 


Hi (o, В,х)+ Hi (a. В, х) 
2 

Н! (о, B,x)- H (a, b,x) 
2i 


e = J (a, B,x)= 
HY(a, B,x)= J (a, B,x)+iY (a, B, x) et Ha, В,х)= J (a, B.x)- iY (a, B, x) e 


ï (aps) 


Et la formule de continuation analytique de la fonction d'onde de Bessel Y < double tildée > : 


x 
Соз(рл)+е 2 
Sin[pz) 


Il vient les deux formules de continuation analytique pour les fonctions d'onde de Hankel : 


J (a, fach e" J (a, B,x) Y (a, В.е") е" (а, В,х)+ 2i Sin[mpz) J (а, Bx) 


= P 
Соз(рл)+е 2 


Ha, B, хе"? )- e" J (a, B,x)+i e "" y (а, В,х)+ 2i Sin(mpr) Sin(pr) 


J (a, p.x) 


LE 

z 2a 

H Na, B, хе" )= JT (a, B,x)-i e™ Y (o, В,х)+ anders 
in[pz 


p 


J (a, B,x) 


P 
Н!(а, B, x) HO (a, В,х) Ha, B,x)- HO (a, В,х) Cos[pz)+ e 2a ( 


Ha, B, xe"*)- e" x pie nmm ; Sin(mpz) Srel HV(a,B,x)+ le, ll 
> p 
| | Hi н? | н -HO E 
Ho, B xe” )= ет p (a,B,x)+ p (a, ,x) ет Р (a, B, x) p (a, B,x) Sin(mpz) Соз(рл)+е É (ut (a. p. x) Ha, B, x) 
Р 2 2 Sin(pr) d ü 
PL nM s 
Sin(pz)H (a, B, xe"* )- -n (a, px Sin((m-1)px)+ Sin(mpz)e 2° - Н (а, В, x)Sin(mpz e" eg 2а 
> 


dub. nes 
Sin(pz)H o, B, xe"*)- itas slo +I)pz)+ Sin[mpz)e Ze |+ itis +e ? | 


Comportement asymptotique des fonctions d'onde de Bessel lorsque a est très grand 


Au premier ordre, lorsque a est très grand, on peut écrire pour la fonction d'onde de Bessel de 
première espèce : 


P р 
232 =; ÈP. Ge? Ea 
2 | e 4 WhinakerM[9,2,,ia 3 2 | ë 4 EE d 
binds 4a 2 " 4a 2 
PEUT T(-p) a x? Г(1+ p) lax? 
lt 
5 i e * ушем A а Ë Я 
ТЕ + NIST 10.167 z=*Â y=2 
2 2 zi P) [04 x? 2 2 


Je dois utiliser le comportement asymptotique de la fonction de Bessel pour un argument grand, 
soil : 


=> œ >> 0 ena) ET — dk Sol 
(22 ee 


Au premier ordre, lorsque a est trés grand, on peut écrire pour la fonction d'onde de Bessel modifié 
de première espèce une formule tout à fait similaire qui conduit à la méme évaluation soit : 


a 1 e ol ape) 


Pour la fonction d'onde de Bessel de deuxième espèce Y tilde, il est possible d'établir que : 


p dei 
em пе E? 
"ER 


4a 


Et comme le comportement asymptotique de la fonction de Bessel Y pour un argument grand est le 
suivant, on arrive à une formule asymptotique au premier ordre : 


1 
SE VT m a x? 2 . fax л(р+1) 
Yiz)se| = | Sin z-=>>=|=> Y x S 
=) 2) d 2 dÉ d 2 | Iz q x? d 2 4 


2р dei 
= a >> 0 EE | j 2 ZÉ xt) 


4a m хуа 2 4 


Avec Іа formule de connexion entre les fonctions d'onde de Bessel J et Y tilde, оп arrive аи même 
résultat : 


2 
, o = 9) 


Еа) ЗЕЕ ИКЕА 


S | de 
| гў E d 2E ET SE rl : d) cal B | Е "E d), 
ч Ë á a. (52) ca (st gÉ ul ol? sta) 


_ ава) É J NA [z СЕ zb и) co sa) 
2 


ES xt Va Sin xa a а 


e B S l I+ zi, (ax all+p)\_ ç B GEN ax rll+p) 
STE | Gi PP К ДЕ “| | м E 


Pour la fonction d'onde de Bessel K, là encore j'utilise la formule de connexion : 


(erar (es A 
p 
I ,(a,B,x)= : М. Е s cn le Esp a) 
eE] a 
4a 
А p 2ле `2) 
edipi вале сорау ду) 5 i | ; I (a. f.x) 
qe "n HE a 
2 4a 2 4 
Qui devient pour a très grand : dÉ 
а к ` sek (n, Вх) = -(1 Cospa))t a. x) = I. (a. p.x) 
T 2 
di 
1а, a Z] — ol Че 
4a (72 je 2 4 
a >> 0 7 
I (a, EE j - Co E s 
40) [ Pa 2 4 
4) (- p) (p+) 
De 4a 2л ах? л(1-р ax? л(р+1 
—$ K (а, В,х)= С 1+C Cos ——— 
SE in(pz) Jo, Bal raw (e d > 1 | (1+ Соѕ(рл)) d 5 1 | 


Comportement asymptotique des fonctions d'onde de Bessel lorsque B tend vers 0 


J'illustre par les deux limites suivantes le comportement asymptotique au premier ordre des 
fonctions d'onde de Bessel J et Y tilde : 


—P 2 
Lim (2) CHEN yx 1 
500 42 E 2 m 2 
2 
aT p+l 
Lim (2) £e B,x)= y 
pec dog Sr KE? P| 2 


On établit ces limites grâce aux deux correspondances des fonctions de Whittaker avec les 
fonctions de Bessel modifiées : 


1 
WhittakerM (0,v,2z) - x Haf v zi, (z) 
WhittakerW (0,v ,2z) =, E (z) 
л 


De même ауес les mêmes formules de correspondance, оп établit les limites concernant les 
fonctions d'onde de Bessel modifiée : 


Lim (2) зела) ys E 


B—0 


B—0 


Lim DEC ,x)= лут 2 


Comportement asymptotique des fonctions d'onde de Bessel en argument x 


La fonction de Whittaker M a pour comportement asymptotique pour de larges arguments z 


2 22 1, -=K 
наке (к,г) DH) is, TH) o 2 E g Jr 
2 
EEN M rats] Comme ü =Z E 8 Eg 
2 4a 
pour u-K* d - "+ô < Arg) Z -ô 
ET Em E 2 E Gees 2 Gs 
De 4 ШОО D iae NUZ E =. | - ee 2 c ) - 
a (2.2) k | 
2 4a 2 4a 
+p or 2 + ax 2 . l+ [È ogla x? kp -i ax B ogla x? 
=Г(1 SCH бү s porcos E 1 x xd Dum | us | ed ) 
Div + pe m p T 2 rí +p) 1+ 2 1 2 
p] P E SE 2202 
4a 2 4a 2 4a 2 4a 


Cela conduit à établir le comportement asymptotique de la fonction d'onde de Bessel J 


р 
232. p EP: 2 
Z) e * taken -iÈ 2 ia J 
a a 


Pour la fonction d'onde de Bessel modifiée I, la formule asymptotique au premier ordre est la 
suivante : 


p 
2 X5 -in t2 2 
Z) e 4 LL d 
4a ` 


S 
= 
EU 


Pour ce qui est de la fonction d'onde de Bessel Y tilde, partons du développement asymptotique au 
premier ordre 
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Avec l'expression suivante de Іа fonction d'onde de Bessel Y tilde : 
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J'obtiens l'expression asymptotique de la fonction d'onde de Bessel Y « tilde » dans le cas de grands 
arguments x : 
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Comportement asymptotique des fonctions d'onde de Bessel lorsque a tend vers 0 


Cette section est dédiée à la vérification du comportement asymptotique lorsque a tend vers 0 des 
fonctions de Whittaker M et W redonnant les fonctions de Bessel usuelles. 


Les fonctions de Whittaker M et W ont pour comportement asymptotique pour de larges paramètres 
x des formules données dans « NIST Handbook of Mathematical Functions, » section 13 
Whittaker's functions (formules 13.21.1 et 13.21.2). Encore faut-il tester ces formules et il s'avère 
qu'elles sont particulièrement justes lorsque le produit x z est constant et un peu moins de manière 
uniforme avec les valeurs de z croissante. Cela nous convient puisque c'est justement la condition 
remplie lorsque a tend vers 0 ! La formule 13.21.2 n'est pas aussi bonne lorsque le paramètre В est 
très grand même si la valeur équivalente de x est de facto également très grande. 
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D'ou soit pour a tend vers 0 : 
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Il est clair qu'en partant de ces seules expressions asymptotiques et en faisant jouer les formules de 
connexion (fonctions d'onde de Bessel et fonctions de Bessel) on arrive aux deux autres formules 
asymptotiques des fonctions d'onde Y tilde et K : 
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Mais voyons si l'on ne peut pas obtenir directement ces résultats à l'aide des fonctions de Whittaker 
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Les limites impliquant les fonctions Gamma se calculent comme suit lorsque le paramètre | è est 
4a 


grand : 


zu 
Ta 2x z ?e? = z=x+iy=T(x+iy)= däs (xiy)" 7а SH 


Haapt, erste ШЕ елү ешге, 


1 
x +y P 


iy 


1 іу+х 2 
РА e Ls 
i KS 


Comme 


p P g 2 2 2 

2 2\2 4a DÉI ре QVE p2\2 — £J 2 E DL B 
[zr [2л ei E e 2 4a zre 4e < Z) В Ë) SEL fre вао “9 É 
(04 (04 


р (p-1)z 
_ B 232 | ШР 2 
К (a, ‚х)=ле sa| В Ке É zx WhittakerW us Ee 
4a) Мах p+1 .B 4a 2 
pr 
2 4a 


Appliquons la formule asymptotique de la fonction de Whittaker K : 
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Comme u= 


Dans la formule de la fonction d'onde de Bessel modifiée K : 
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= K,(a,B,x)= "rg ie U B x)+iY,(ip 2) - "rie Hp 2 


Cela donne donc : 


Ce qui compte tenu de la formule de connexion redonne bien la fonction de Bessel K : 


= K,(a,B,x}= Erg ié Hs s) = Re(K ( x))= K (8 x) 


Comportement asymptotique des fonctions d'onde de Bessel lorsque B très grand 


La formule de NIST 13.21.1 s'applique bien lorsque B est également vraiment très grand, ce qui 
donne l'évaluation précédente : 
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Les deux formules de NIST 13.21.1, 13.21.2 ne sont pas totalement satisfaisantes notamment pour 
les fonctions de Whittaker W lorsque l'argument z est constant et B grand. Partons de l'expression 
asymptotique de la fonction de Whittaker M : 
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Cela nous donne donc les expressions successives suivantes : 
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Је пе parviens раг à simplifier suffisamment cette expression ауес les paramètres correspondant 
aux fonctions d'onde de Bessel. 


En revanche j'ai constaté que les deux formules de NIST 13.21.6 et 13.21.7 sont bien meilleures 
d'application lorsque В est modérément grand. D'autres part ces deux formules ont le mérite d'être 
relativement simple pour les deux fonctions de Whittaker. Elles s'appliquent tant sur les paramètres 
originaux des fonctions de Whittaker que sur les paramètres liés aux fonctions d'onde de Bessel. 
Les voici : 
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Une forme équivalente est la suivante : 


èla) _ Vz+2° + Logl/z + free WhittakerM (- K, UA z) z V2 rens 1) <, (4c (z)) 
2 2 


SE T WhittakerW(- Ak :)= Z> |®\© g, (акё (2) 
Z 


Avec les paramètres utilisés pour les fonctions d'onde de Bessel, cela donne pour la fonction de 
Whittaker M : 
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Et l'approximation asymptotique de la fonction d'onde de Bessel est la suivante : 
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Pour ce qui est de la fonction de Whittaker W après de multiples calculs pas si complexes que 
cela : 
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Ге résultat s'avère finalement plutôt simple : 
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Pour trouver l'expression asymptotique de itae i£, ia i j'ai consulté un article de 
a 


1957 des deux auteurs A.Erdelyi et C.A.Swanson < Asymptotic Forms of Whittaker s Confluent 
Hypergeometric Functions », dans lequel il donne la forme asymptotique de l'expression de la 
fonction de Whittaker M (formules 5.5, 5.7 et 10.1 ) 
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Avec les paramétres utilisés pour les fonctions d'onde de Bessel modifiée, cela donne pour la 
fonction de Whittaker M : 
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Et l'approximation asymptotique de la fonction d'onde de Bessel modifiée est la suivante : 
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Pour la fonction de Whittaker W a lieu le développement asymptotique suivant au premier ordre 


(formule 10.2 du même article) : 
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Soit en appliquant sur l'expression de la fonction d'onde de Bessel modifiée K, et en intégrant la 
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formule approchée de la fonction Gamma suivante d 
2 4a 


p (p-1)z 
-E( BN2 | "a 2 
К (a, ,x)=m e S E Re 2 Whittaker W ;, Р ija x? 
4a Jax r p+1 PU 4a 2 
2 4a 
(m rei 17 B2Vaa 2 
e * 2e E j A i р e ) 
q Г =, 
1 ¿(z pi DH jP" Е 2 V 4a 
Re | 2 PAN БУЗ 
2 4a 


K (Bjr er CG) 


1 2 a 2 S 

(2p-1) 109 рт _ B 2 үз (uf 
4 E x2e 4 a 2 | 2 

2 4 2 

"y queo pe 

2 4a 


e X ———————— 

1 c (z 1 8° P 
К (a, p, R Sg. 
„(a B x)= л FE \ 2 (2) (z ) e 2 


Les limites suivantes représentent le développement à l'ordre 0 de В (soit le premier ordre) : 
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Le résultat final au premier ordre est également assez simple : 


Dans ce méme article de 1957 < Asymptotic Forms of Whittaker S 


Confluent Hypergeometric 


Functions », les deux auteurs A.Erdelyi et C.A.Swanson donnent les formes asymptotiques des 


fonctions de Whittaker M er W (formules 5.5, 5.7 et 10.3 et 10.5) : 
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L'application de ces formules à la fonction d'onde de Bessel modifiée I donne : 
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L'application de ces formules à la fonction d'onde de Bessel modifiée K donne : 
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Comportement asymptotique des fonctions d'onde de Bessel lorsque p très grand 


Les fonctions de Whittaker M et W ont pour comportement asymptotique pour de larges paramètres 
p, les expressions données dans NIST < Confluent Hypergeometric Functions, Whittaker 
Functions », formules 13.20.1 et 13.20.2. 
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Mais j'ai vérifié leur adéquation avec un test numérique et ce n'est vraiment pas satisfaisant. 
Utilisons les tout de même pour obtenir une formule asymptotique grossière lorsque le paramètre p 
est très grand. 
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Pour la fonction d'onde de Bessel Y tilde, il vient : 
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Pour avoir une meilleure estimation, il faut consulter l'ouvrage de 1974 de F.J.Olver et W.Rheinbolt 
« Asymptotics and Special Functions », chapitre 11 sur les fonctions de Whittaker , oü sont données 
deux formules au premier ordre en page 261, comme suit : 
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Ces deux expressions assez complexes, peuvent se décliner tout de même avec les paramètres 
utilisées pour les fonctions d'onde de Bessel. Commençons par la fonction M de Whittaker : 
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Cette expression peut encore étre simplifiée car l'expression de valeur complexe sous l'exposant est 
en réalité de module 1. En effet : 
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D'une part on retrouve bien la formule ultra-simplifiée donnée par NIST, grâce à la limite 
suivante : 
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Et d'autre part dans la fonction d'onde de Bessel J, il vient : 


P xP 2P 2e 2 2 
(a, B,x)= —P p 


2"T(1+ p) er ete в) о а т) 


Avec la fonction d'onde de Bessel modifiée I, il vient (il suffit de remplacer dans l'expression de J, 
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L'expression de la fonction de Whittaker W avec les paramètres utilisés dans la fonction d'onde de 
Bessel Y tilde donne : 
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Les fonctions Gamma peuvent également se développer asymptotiquement comme suit : 
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Comme on peut appliquer toute une série de limites lorsque p est très grand, on arrive au même 
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résultat de la formule simplifiée de NIST : F (a, p.x)* Lg 


En effet : 
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Pour Іа fonction d'onde de Bessel modifiée К : 
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On retrouve bien l'expression asymptotique simplifiée puisque : 
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Représentation intégrale des fonctions d'onde de Bessel J et I 


À partir de la représentation intégrale de la fonction de Whittaker M, il est facile de construire une 
représentation intégrale des deux fonctions d'onde de Bessel J et I. 


En effet la représentation intégrale de la fonction de Whittaker M sur la ligne réelle s'écrit : 
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Avec les paramètres des fonctions d'onde de Bessel J, il vient : 
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Avec Іа fonction d'onde de Bessel modifiée I, on arrive à une résultat aussi simple à l'aide de la 
formule de connexion. 
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Les conditions pour lesquels cette intégrale est définie sont Pour Dr E pue asc Cette 
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E 
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‚ү 
conditions assure l'intégrabilité du terme (-! ) ` sur [-1,1]. 


Ce n'est pas aussi simple avec la fonction de Whittaker W et donc pour déterminer une 
représentation intégrale valide de la fonction d'onde Y tilde. 


Si je prends la représentation intégrale de la fonction de Whittaker W sur la ligne réelle : 
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Pour ek > Re(x) L4rg(z) < - 


que l'on peut également écrire par le changement de variable d'intégration t en 1/t 
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Pour Без > Re(x) L4rg(z) < 2 


Alors elle n'est pas applicable car l'argument de z devient justement égale à n/2. Et pour cause, car 
dans ce cas l'intégrande devient une fonction oscillante dont l'amplitude est croissante. Je me suis 
amusé à calculer explicitement la représentation intégrale de la fonction Y tilde, méme si l'intégrale 
n'est pas réellement défini. On obtiendrait formellement : 
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се qui confirme bien que l'intégrande est ипе fonction infiniment оѕсШапіе d'amplitude croissante 
sur le domaine d'intégration. 


Avec les fonctions d'onde de Bessel J et I < double tildée », il vient : 
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Autre représentation intégrale avec un noyau de fonction de Kelvin 


On peut également utiliser la représentation intégrale suivante avec un noyau de fonction de 
Bessel : 
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Avec les paramètres des fonctions d'onde de Bessel J et I, il vient : 
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Formules de connexion des fonctions d'onde de Bessel J et Y tilde en ordre entier 


Avant d'utiliser les formules de connexion déjà établi, et d'opérer un passage à la limite lorsque 
l'ordre devient entier, voyons quelques expressions impliquant les fonctions Gamma. On sait que les 
expressions suivantes des fonctions Gamma avec des ordres p entiers s'écrivent comme suit: 
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Cela va permettre de définir plus précisément le passage à la limite à opérer lorsque les ordres sont 
entiers afin de définir la fonction de deuxième espèce : 
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Pour les formules limites en ordre entier, le développement à l'ordre 0 des expressions ne suffit pas, 
il faut pousser au moins à l'ordre 1 : 
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Ce qui donne finalement sur les fonctions d'onde de Bessel de deuxième espèce : 


E 


p =0= Y (a, В,х)= Lim - : — S кнн x)-J (@,B,x) 
B 2 4a 


yn Sin(vr) а 


v>p Sin 


e 
1 k=p-1 4a? É 
p pair = Y (a, B,x)= Lim TA (56 | Cos(vx )+e 
J 


т я ai 
Cos(vr e" +e ** 
1 LE 
4a 4a 
эй D E eon e =€ 


efist- deal"? kelly: tn tan 


EN 
Cos(vz)+ e a ee [nest Jn 
m m ue p 2 4a 
p impair > Y (z, B, x) = Lim 


vo» Sin(vr) NE 


Ties spot: ege 


On peut encore simplifier ces formules en négligeant dans le passage à la limite le terme 
comportant le Cosinus, soit : 
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Les formules de passage à la limite illustre parfaitement l'application du cas limite a->0, où l'on 
retrouve les formules des classiques fonctions de Bessel, à la condition que : 
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Lim Log . +Re w EE =0 
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Ce que l'on peut vérifier numériquement, par exemple, grâce à Mathematica. Le comportement 
limite de l'expression est le suivant : 
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ce qui valide le passage à la limite a->0 vers les fonctions de Bessel traditionnelle. 


Dans le cas où a est trés grand, soit Bun trés petit, la formule de connexion limite pour les ordres 
entiers prend la forme, en utilisant les comportements asymptotique suivant : 
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Il vient les formules limites suivantes de connexion : 
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On retrouve bien la relation de connexion limite pour les fonctions de Bessel J et Y asymptote, étant 
donné que le terme (v-p) l'emporte sur le terme en logarithme des paramètres a et P. 


Démonstration des deux conjectures sur les formules de connexion limite pour les ordres 0 et 1 


Nous sommes maintenant dans la capacité de démontrer que: 
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Prenons les deux formules limites de connexion établies : 
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Et prenons d'autre par t les formules de décalage des fonction Y et Y tilde obtenues pour les ordres 
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Partant des formules limites, il vient : 
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En utilisant les deux formules de connexion, il vient immédiatement le résultat à démontrer, soit : 
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Ces formules limites de connexion nous permettent d'imaginer des nouvelles fonctions de deuxiéme 
espèce qui possèdent des formules limites de connexion encore plus simple : 
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Soit en introduisant les deux fonctions suivantes : 
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Les nouvelles fonctions introduites ici sont identiques aux fonctions de Neumann dans le cas a=0 : 
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Dans le cas a très grand, il vient : 
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Simplification де Іа formule de connexion des fonctions d'onde de Bessel J et Y еп ordre entier 


On a vu que les formules de connexion des fonctions d'onde de Bessel en ordre entier se présentent 
comme un processus limite : 
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Mais l'on peut encore simplifier sous la forme unique suivante : 
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Une fois déterminé le décalage entre les fonctions d'onde de Bessel Y et Y tilde, nous serons à méme 
de donner la formule de connexion pour les fonctions d'onde de Bessel Y telles que construites par 
Moon et Spencer. 


Formules de connexion en ordre entier pour les fonctions d'onde de Bessel modifiées I et K 


Si l'on revient à la formule de connexion des fonctions d'onde de Bessel modifiées : 
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Transformation de l'équation de propagation de Bessel 
La transformation la plus simple consiste en un changement d'échelle de l'argument : 
(6) Equation de propagation de Bessel Z'(z)«— 124 )+ {= z? +q- 5 ke =0 
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Mais comme en même temps la solution de l'équation peut s'écrire : 
V(yz,z)= AJ, „(x,1,z)+B,J_„(z,1,z) 


On en déduit une loi de changement d'échelle pour les fonctions d'onde de Bessel : 
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Par la transformation z->1⁄z l'équation de propagation devient : 
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On observe que l'équation possede maintenant un póle d'ordre 6. 
Lien avec les fonctions paraboliques cylindriques 


La transformation suivante est légèrement plus complexe : 
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Par ailleurs il s'avère que cette même équation résultante est également de la famille des équations 
de Weber dont la solution est la fonction parabolique cylindrique D : 


V "(z)+(a?z? +в V (z)= 0 


Posons z=w = V"(w)- 2 (o? 2w? * BY (w)= 0 


=> V'(w)+ 2(0° 2 ур? +e V (w)= ZOnE ryo- 0 


4 
1 i 1 (1+2) 
> Va =--> 2 =+=> A = = 
4 2a (i-iVa 2a 
.02 .о2 
pp sn APE е n а 
2 2a 2 2a 
= - A2 oi 
Choississons 2 l D 1.294 gs у= 118 — 1-V= 1 (8 
24a À (1- i) 2 20 2 2q 


Solution V.(w)= D (w) V, ў =D,,(iw)= 


tei ka) v, ()-D , „(еъ дема) 


“ Lip” HP 
2a 2 2a 


2 


On peut donc établir un lien formel entre des combinaisons linéaires des deux types de solutions : 


кө (a fx) =D ftne) 


2 2a 


OEA В) V;(z)=D | jp (iia) 


2 2a 


De méme pour des combinaisons linéaires des deux types de solutions suivantes : 


: B' 1. uos 
WhittakerM ‚ех 
4ia 4 
Jx 
2 
WhittakerW В ETT 
4ia 4 


dr 


V (x)= BAU V (x) = 
b =i Je: Ba)" 


V, (x ) = 
et donc finalement pour des combinaisons linéaires des deux types de solutions suivantes : 


2 
Winakert L Э ,ia x 
V. (x)= 4ia - 
| A e [E 
iieri L Lues V(z)- D MED 
ia 


dr 


bic a) 


En consultant la base de connaissance NIST sur les fonctions spéciales mathématiques, à la section 
« relation avec d'autres fonctions spéciales », les liens sont établis explicitement : 


1 
WhittakerM LE e Ld =2* T i32 [Ett tz) ez T Lus E D i a,z |+ D he 
2 42 2 4)\л 2 4л 2 2 
1 
WhittakerM E = 2? LÉ ET z) U(a,z)}= 2? LÉ | H | a, =z H Í a,z 
2 42 2 4)\үл 2 4)\л 2 2 
š 1 1 1 2 za T 1 
WhittakerW = GECE =:22 JzU(a, z) - 22 VzD ce Z 


On peut donc explicitement établir un lien entre ces fonctions : 


2 2 Fe 2 2 
Love ER => | а= L dp li l=-1 В L 2==2 B 
2 4ia 2a 2 2 2а 2 4ia 2 4ia 
1 ,3_3a+ip? 1 2 1 a+iB? 
2 4 4a 2 4 4a 
11 las (1 1 1 
WhittakerM ED y |= 22 ER н (2 D a,y |+ D| -—-a,-y 
42 2 4)\л 2 2 


1 
iax = y y =2ia x — y =J2ia x= (1+iWax= Jy =( ech 


WhittakerM сере 3 " "P Së an | CARE Tn RET E $ d Ee Nas) 
4 4 2a 2a 


ia 4а Ут 
i - iter LZ s) Puer E Lias) 
Comme J (a, В,х)= В Ss = Ï i 
3 d xia) e 
2 e Jo | a +iB? eMe en оа) 
E (a B Se 4a 2a 2a 
A Dor 


1 11 122 (1. p 1 l 
WhittakerM a,—,—y° |= 22 T a+ Ур a,-y - D a,x 
2 42 2 4 )\ x 2 2 


2 EN “22 А n WER oi 
ЕЕ er = Ehf E (ed 799 eis] 


ia 4a Jn 2a a 


2 
s ; ишен, ла 3 ne d pinaten P. ie d 
Comme J (a, В,х)= B Ge EE Е 


2 d x(ia Ул x(ia)a 
2 gean | CHE , Маз d arif E, 
io ,x) á : 5 - 


Les graphes comparés des deux fonctions ainsi calculées, confirment l'égalité : 


= BesselJWaveRecurrentSeriesFunction[a, 8, L x) 
2 


g 
zaia уй ik =a 


Pas x Va P au We * va) 
\ 2a 


Re 


n2" H 


a = 3.5; 
В= 1; 


— BesselJWaveRecurrentSeriesFi unction(a, B.- Е х) 


E (3 a-i 82 
4 


Lm o genie yo (14) x Va) 


ET 


EVE 


SC x VB 


Dans le sens inverse de la transformation de l'équation, on obtient bien le même résultat : 


Equation ` Z''(z)+ (p° + a?z")Z(z) =0 

Posons Z(2)=Vz U(z) 

U(z) U'(z) „© HG 
2Vz We Dus. 4 


Z"()+ (022° + B*)Z(2)=0=> VzU"(2)+ == +(a2z?+ в? Wzu(z)- 0 
5 (e: 


? +B? - Tyo- 0 


U (z) solution de l'équation de propagation de Bessel 


Ce IDE E Vo- 0 


Z'(z) = 20 (2) + Z'(2)=VzU"(z2)+ 


eU —— 


en posant U"(z) + 


> p= Solution Z(z)= A.J (a, В,2)+В,Ј (о, B,z) 


2 2 


En inversant le signe de l'un des paramètres, il vient : 


Equation Z''(z)+ (s? - a? z? zc) -0 

Posons Z(z)- VzU (z) 

U(z) 
24 z 


Z"(z)+(a2z* + B*)Z(z)=0=> Jzu' (2 


UG), Ui Об) 

FE 2Vz 4zJz 

U'(z) ` U (z) 23 2 = 

f cre z? + B? WzU(z) = 0 
2 2 __ 1 ú 

Je — OT zo 0 


U (z) solution de l'équation de propagation de Bessel 


Z(z)y2JzU' (z—-— Z'(z)5JzU"'(z) 


EUM TEM 


en posant U" (z) + —— — Ju +В? = Z Vo- 0 


> p=> B—iB Solution Z(z)= A.J ,(G,iB,z)+ B.J ,(a,iP,z) 
2 25 


Га dernière transformation que nous étudierons est Іа suivante, elle confirme Іа transformation 
précédente : 


Equation Z''(z)+ (p° + a?z?)Z(z) =0 

Posons 2(2) = 2 О(А 2) 

Z'(z) »z'AU(Az)*yz""U(Az) 2"(2) = Si P#U"(22)+2y Az" U'(Az)+y(y—1)z”2U(A z) 
Z'(2)+(a?z? + B*)Z(z2) =0=> z" А (А ele Ae â z! U'(Az)+y(y —)z'2U(Az)+[a222 + 8° ku z) = 0 


—U"(A ae ZU de аг yp 10 Zt ya г)=0 


2 — 
w-Az —— — pare D pee 
w w 


Des formules liant les fonctions de Whittaker aux fonctions paraboliques cylindrigues, on peut tirer 
un autre expression de la fonctions d'onde de Bessel J : 


Pour les fonctions d'onde de Bessel Y tilde, I et K, il vient ; 


y = “pl 2, ,; 8° Ll Boss 
"isle dept LE NM — г) 
(a B,x) rh. GÉIE | EE 
4 4a 2: `2 
dr „B D 
E. Be * 2 же? | 2 | 
I (a, B,x)= R Dc 14 ia. 
GEES - 
4 4a 
Tagh: SE o 
6, ,x) SC Re Al iim (is) 
4a 


4 4 
ы p D iZ 2 
Pre © оеро 1 B 
K ,(a,B,x)= e S d +i (reis 
AE d a ET 
4 4 


Z 


2 2 2 2 
40-4219 jee 2 | ou Z(z)= AJ É Lan É | 
E ON. 2 +562 EX 
2 2 2 2 


1 Se a z? a z? 
q=0 p=0=>Z"(z)+=Z'(z)+a?z?*Z(z)=0= Z(z)= A.J, + B.Y, 
z 


i-o ze) rep (az bi 


gees ech Lotte Zu 


= Z(z)= AJ,(qz)+ B.J_ (qz) ou Z(z)=AJ,(qz)+ B.Y, (a z) 


Quelques formes équivalentes de l'équation de Bessel 

Voici quelques transformations de l'équation de Bessel permettant d'exprimer les solutions 
d'équations différentielles apparentées : 

1 2 p 1 p 
()Z''(z)+—Z'(z)+| B -= Z(2)=0 QZ'(2-—Z'(G)|1-— |Z(z)= 0 

Z Z Z Z 
Depuis (1) en posant Z(z) -V(Bz) z-BzZ'(z) = B'V"(Bz),Z'(z) = BV'(B z) 

1 2 

DS BV"(B2)+_BV'(B ole Con z)=0 
peas i ria e 

B z B? z? 
Z(z)=A,J (B z)+B,J_,(B z) 
/6/B 5, is zq- > „VG B) 

Z 


a+l 
Z 


yu z)20— V(z)- AJ,G)-* BJ ,(G) 


Depuis (2) en posant Z(z) = 


Hi E s VEB), 
Biz E 
DEEP a EP aan EA, (ZEA a EP), Ë Con. 


B? z* B P 2% B z” zet! 


V(z/ p) 


Z''(z) = ala +1)=>5- „ar 


25 1) 


V"(z/ V'(z/ V(z/ 
= (z/ B)— neq aye 


= V''(z')— — t EE ye 0 


2a -1 
Pour l'équation V" (z)— X= 


Wé (rss ya 
=> solution V (z)= z^Z(B P V(z)- z (A.J, (8 z)+B,J_, (8 2) 
ИЕТ 

2| = p" (z ` Е да) СЕРТИ) 


2 


[= p= 
P 


Continuons par une autre transformation de l'équation de Bessel : 


Depuis (2) en posant Z(z) = GPO z'=(z/ p)” => z = B z” 
pray ZIE aED 
B e pe 
пу . Ву?) gä E АІ ЕМА АС BI i ooi) Ee BY) 
po P yz B'yz z 

220/7-) V"((z/ py”) А "m i 2а) zur V'((z/ BX”) t i V(z/By") , 

By? z^ В'?у д9 z^? 

zur oe BY") o V BY"), P VEAD _ 0 

By P P z? 2° 

n E C s ns “+ (D) Jup -p | VED 0 
B ''y z В yz Z 2 


> —————— 
py z^ ET 2(/y-1) 


By V'(z) a =p? V(z)p^"y? 
=> V"(z')+(1/ 2a) Ber 1+ 2 2(/y-1 =0 
= (ez) | re | 


' (> E p^ 2 
=> V"(z')+(1-2a y) Уе [pe Pepe ë jee 
Z 


z 


2?(/7-) (meur 2a) BY Hi ((z/ B?) , Jup FRERE )- 0 
z? 


Comme Z(z)= A.J, (z)+ B.J_, ( z) 
= V(z')= (8 z" Y (4J „(Bz")+B.J_,( B z" ))x SE (4.7, (в z" )+ в.2 „(В z") 


E 
Solution de l'équation V'(z)+(1-2q y) ——— — dÉ Aa LE ZP E 
z 


= V(2)=z"(4,J,( Bz')+B,J_,( Bz')) 


Avec ay >a 
A ИРЕ 
Solution de l'équation V''(z)+ (l 20) ——— o. Je EE + e yo =0 
z 


= V(2)=z'(4,.J,( P z")+B.J_,( B z" )) 


Résultat 
252 
(4) 2"0)- — "ele 2g. PL bei 


Z(z)- z'(AJ,(Bz')- BY (Bz) si peN 
z()-z"(42,(Bz')-B.J (8Bz') sipeN 


1 2; 2 2 2 2 1 1 
En posant a =—;y=2; p = Ма; = È => e È! В = 
P 2 y y p 1 p B 


(4) > Z''(z) + X'a°z°Z(2) = 0 


= Z(z) = Alle н 8.7 |2 az ) 


En partant de l'équation de Bessel : 


Oz"o«lzo Ë E ре =0 
2а -1 
Z(2)=z"(4,.J (Bz')+B.Y (B z) si peN 
Z(2)=z*(A.J (Bz')+B.J_.(Bz')) si p eN 
EE 1 


1 1 
En posant œ == ; =—;y'=2y-2 у= ;p=== BY = a 
p SPY dd pa 27 en By 


L'équation (4) Z''(z)— 


2 25512 
ze» Byz”? ET ku — 0 
Z 


a pour solution | 


a"? 2g? 


Y B y'+2 


(4) > Z'(z)+a'z"Z(z)=0 a pour solution 


d. 1/2 y'+2 1/2 у'+2 
H A E м 
s WA SW Д 


y'+2 


En ôtant les primes cela donne bien 
(5) Z''(z)+aœ z”Z(z)= 0 


1 1/2 у+2 12 у+2 
> 2 L 2 у+2 
a pour solution o=" 4 | | Z 3 z? E | | EE ) 2d 
m + wg 


y+2 


eN 


y+2 


y+2 


1 1/2 у+2 1/2 y+2 
Í Зат A VY pe | 
ouZ(z-z!| AJ, | Cz? |+ву, | 32 || si——eN 
= 7+2 == yp y+2 


y+2 


Posons l'équation ayant une forme plus générale et montrons que les solutions sont construites à 
l'aide de fonctions de Bessel : 


(8 z'U"(z)«Qa-2v +1)zU'(2)+(B er +a(a-2Bv)U(z)=0 Роѕопѕ U(z) = z"'" Z(y z") 
ауес Z(£) solution de Ç*Z'(Ç)+6 Z (C) (P? -v?)Z(6) =0 


Vp eg VB 
menn ne ЫЕ 
y У y 


a—pv41-f 


re) (ET CEE 
e S CE da EN 
A5 NO ^ rene 


а-Ву-28 
gig 228009 - Bv = “л B 


Con Чү deen 
te SE le - Bg dE PON d с jJ 


B 
«i| (а Bv Xa - Ву E 


; e 
Z"(f)5— 2 
G Z"(6) e 


a-Bv+1-2B 
„apeti 22) B 
By y 


1/8 
+ (2a -2pPv мя) U 
y 


1/8 
S vam? |, 
үү Y 
ES roro RE 


1/B VB 2/B 1/8 
3 — + (2a – 2 Ву + (D 2l I B (8) | 
y y y y 


«a - Bv) «(C г-у (Е a] l 


B 


F= 
= =X 
Ke 
— 
= 
мез = dea 
+ 
A = =>: 
`< |S 
SE 
2 
= 
S 
LAN 
TE 
SR 
` ^L 
ES 
му 


PrE EOE- OLE) 


e EECHER 
= IEN e) + Qa -2fv + (DEN +(a(a-2Bv)+6? ewe -0 


Comme, il vient la solution de l'équation différente cherchée. 
e? = у? 728 
=> (9) U" (z)+ Qa -2Bv +1) zU'(z)+(a(a -2Bv)+ B^y? z^ \U(z)=0 
c.q. f d 


AJ, (y zÊ) + B.J B 
Solution EE JO z )+B.J (у= ) | 


ou A.J, (y 2%) + BY, (y z^) 


Formule de récurrence des fonctions de Whittaker М et W 


Voici quelques formules de récurrence des fonctions de Whittaker. Sur le site Web Mathematica des 
fonctions mathématiques, il est donné la formule suivante : 


16(u -1)uQu -1((2u -1(2u -3)-2zv) . 
SE 37 xm = Yu- 1-2) WhittakerM (v, u —1,2)+ 


2 
16(2u 1) (u —Un WhittakerM (v, u —2,z) 
Qu-1-2v2u-1-2v) 
2u+1+2V\2u+1-2v) 
1би(и+1)(2и +1) 
(Qu «Xu -1)-2zv) WhittakerM 
z(2u _ 1)(2u + 1) 


WhittakerM (v, u, z) = 


WhittakerM (v, и—1, 2) ( WhittakerM (v, u+ 1,z) = 


(v.u.z) 


De méme Pour N.Y.Vilenkhin « Special Functions and the Theory of Group Representations », 
Chapitre 8 Whittaker Functions, page 416 et suivantes , il est donné : 


. 2 
(2u-1) _ a nz) + = 1) v ateq, u,z)=2u(2u-1WhittakerM(v, u -1,z) 
2 2 


2 

: 2 та -v 

(2u+ p OWhittakerM (v, u, z) TD v WhittakerM(v, u, z)= СЕ WhittakerM (v, u +1,z) 
д 22 2(2u+1)(u+1) 


ôWhittakerM (v, uz) + É — 1) 4 rate H, z) = 2uWhittakerM (v, и—1, 2) 


ER 2z (2u = 1) 
D | 
OWhittakerM (v,u,2) + | Qu us 1) Y. WhittakerM (v, Ш,2 ) E Cl i 2 
à 22 Du mu 
‚ i i 
S ƏWhittakerM (v, ш) (Du iu 1) SÉ WhittakerM (v , u, z) = E WhittakerM (v, u+ 1,z) 
22 2z(2u +1) EE 
(2u+1-2v)2u+1+2v) WhittakerM(v, u +1,z)= 


WhittakerM(v, u -1, 
л. 16(2u +1) u(u +1) 


(Qu - Qu +1) ey len ren 
z(2u-1)2u+1) 


(v.u, z) 


Les deux formules de récurrence coïncident donc. 


Sur le site Web Mathematica des fonctions mathématiques, il est donné la formule suivante de 
récurrence des fonctions de Whittaker W et de l'expression de la fonction d'onde de Bessel : 


4и и Län 3) 2zv) WhittakerW 
z(2u-3{2u-1-2v) 


+ (2u-1)2u-3+2v) WhittakerW (v, u =2,z) 
(2u-3)2u-1-2v) 


WhittakerW (v, u +1, z)- er 
H Kescht 


(v.u.z) 


(v, u- 1,2)+ 
WhittakerW (v, и, z) = 


(2u+1-2v) 
а (2u +1) 
= «uu ы 2u ш 1)- 2zv) WhittakerW 
z(2u -12u +1) 


WhittakerW (v , u —1,z)— 


De méme pour N.Y.Vilenkhin < Special Functions and the Theory of Group Representations >, 
Chapitre 8 Whittaker Functions, page 416 et suivantes, il est donné : 


. 2 
(2u-1) Orate e uz) " = 1) ООО u,z)= 4; +4- T rer u-—1,z) 
Z Z 


Qu+1) OWhittakerW(v, u, z) = +1) Y 


WhittakerW (v, u,z) = Ë - u- 2 WhittakerW (v, u +1,z) 
02 22 2 


(2и ы an s 1)- 22 WhittakerW 


4 = 
" 224-122 +1) бл) 
_ Он+1-2) WhittakerW (v, u +1, z)— (Qu-1+2v) WhittakerW (v, и —1,z) 
(24 +1) (2u -1) 


Là encore , les deux formules de récurrence coincident donc, ce qui est heureux ! 


Estimations_asymptotiques_des fonctions de Whittaker de l'article d'Erdeliy et Swanson, 1957 


« Asymptotic Forms of Whittaker s Confluent Hypergeometric Functions > 


L'article en question explore le développement asymptotique des expressions suivantes mettant en 
jeu les fonctions de Whittaker : WhittakerM (xK, 11 ,4x z) WhittakerW (к, u,4x z) WhittakerW (=x, u,-4x z) 


Je me contenterais des développements asymptotiques mettant en jeu uniquement les fonctions de 
Bessel, car l'article présente également ceux utilisant des fonctions d'Airy. Pour cela sont définis 
les fonctions ¿(z) du paramètre z suivantes : 


zelo] o jaa a ER ep] de 210 RED 
eo 
z«0 Am()- ¿(=+ E ARE ces = z<0 goja cines does lt 


» 2l z 


Les expressions asymptotiques sont les suivantes, données directement comme résultats sans calcul, 
ni démonstration (en suis-je pas capable ?!) : 


WhittakerM (xK, u ,4x 2) = Dei E 1) | SCH (4c (z)) 4e (x) < 


K 2 


T 
5 —o<z<a<l 


1 
Whittaker (к, HA z) x nii 19 


{бїї(л(к al. (Act (z))- Соз(лт(к - ш), (4x6 (а) n 


—-oc«z&a«l 


T 
Whittaker W (— k,u,-4K z) = Jak aM SE) (дег (2) wë 


ee) 7 


Premier cas : z < 0 


La paramétrisation proposée de l'argument z est la suivante : 


AS . . 2 . . LT 
z - e" Sinh (c) Ç [e Sinh? (c )- e ? EE Kee =e 2 r Sinh[2o)+ а 
d + іл = H . E 
E men . Cotank(o) - conils) ou C'(2)=-22 ES =-+= 114 - : € Cotanh(o) 
d do \do 2 2 2V z  2V Sin (с) 2 


Cela entraíne l'évaluation asymptotique suivante de la fonction de Whittaker M : 


j K К2) х Tu 1) 42) Lab e Ebene À 
Vitiis aes] ТСЕ! E (acte) œ Ale LEE celo npe 


K 


WhittakerM (c, u,4x e" Sin (c))« fan rQu He" e í^ rose ) ow: {Sinh( 2с)+ 2с}- era 4) 


WhittakerM (c, 11,4x e" Sinh (с ))= fan rQu He” 4 ` fTanh(o ) ouis ÍSinh[ 2с)+ 20}- Zr 2 


Pour la valeur précise de l'estimation asymptotique nous retenons : 


WhittakerM (re, и,Ак e*" Sin (o ))= E r(2u +1)к” E ë < Tante ) ouis {Sinhl 20)+ 20}+e Si ul + a 


Et l'évaluation suivante de la fonction de Whittaker W : 


Whittaker W (к, HA z) x Max Zen — Ç ) Ísin(z (x = u)J 2u (4c (z )- Соѕ(л(к = vllt, NUS (2)) 


л) s e ae) Н 4) 


ele nd st del 


Sin(r(x — ден а d С al 
- Colle dich) (^ x d 


=> WhittakerW (к, u,4K z) x V2x"e* 


FG 
=> WhittakerW (x, H. 4ке" Sinh ( c)}= 2k*e SS ч кайи аш z)+ der zez) 
=> WhittakerW (x, 11,4x e" Sink’ ( c) 2к*е gt 4 + [Tank sexa z)- de | 


> Whittaker Ws, ‚Аке? Sinh? ( c) "e "e * "regel: {Sinh(20 )+ ele 56-2 
Pour la valeur précise de l'estimation asymptotique nous retenons : 
Whittaker W к, ш.4к e" Sinh? (o))= 2x"e "e * fT ano os (silo): 20}+ ex = a 


Etant donné que les deux expressions suivantes sont également valables avec une valeur complexe 
de x d'argument -л/2 : 


WhittakerM (s, HEI? еі (c))« Em: + leie? Тапс EE 20}+ era + a 
T 


Whittaker W (x, Odre" Sinh (с))= 2к*е se * + ае Z K{Sinh(20 )+20}+e? a) 


En posant к=-іт, il vient : 


T 


SE 
xT(2u+1)e? ү 3) ‚у Beca s (Sint(20)+20)-n (+L) 
T 


1 
> WhittakerM (ie, ir (2Sinh(o) = r(2u+1)e? ak d E SEN elle Al 


T 


pt D s... M ise eaeleis[e-1] 


WhittakerM- iT, H, ir(2Sinh(o ND г(2и+1)(-+ ef? "y FM Los совин о) чечан 


AT 


m 1 
sie Zell D Tanhlo хато | Sinh( 20) +20 len + a 


N 2e об) JTanh(o c с )Cos cl ÍSinh( 20) +26} HE | E sien] 


WhittakerW(- it, H, ir(2Sinh(o y Je a orl- (0) е" E "tie" Sweet, 


AT 


=> WhittakerW ic, uic (2Sinh(o y = e? elrlt-1os(z)-sim(2o)2o) танно) 


D'ou une estimation asymptotique des deux fonctions de Whittaker dans ce cas : 


REN 
WhittakerM(-it, uic (2Sinh(o)Y )= r(2u+1)e? | es [RH euch A 
T 


AT 


WhittakerW ic, uit (2Sinh(o ) L е? зол Log(r}-Sinh(20)-20} (rank) 


L'une de ces estimation asymptotiques est effectivement donnée dans l'ouvrage d'Erdélyi et 
Bateman, sur les fonctions harmoniques du paraboloïde de révolution, « Higher Transcendental 
Functions vol II, 8.7 Parabolic Functions and functions of the Paraboloid of Revolution > formule 
1l page 127, se référant à un article d'Arthur Erdelyi de 1937 dans un bulletin de l'Académie des 
Sciences de Vienne (en allemand), mais introuvable à ce jour sur internet par exemple. Toutefois 
l'article de 1957 d'Erdélyi et Swanson donne un résultat équivalent. 


Deuxième cas : z trés proche de 0 comme 0(к23) 


On évalue directement la fonction ¿(z) par sa valeur approchée : 


с) * deal + ArcSin| z ) s BS 


et s'er => Sen 


i 20 
WhittakerM Le, u,4x z) = j TE 
DER vi? 


> 


Y, (4x6(2))= Y, [ac =) 


EM g 
WhittakerW (к, u,4k z)« 2/24 nk ?е*Мг : 
а Е 


Remarque : c'est exactement се type de cas qui permet de justifier de Іа limite des fonctions d'onde 
de Bessel modifiées vers les fonctions de Bessel modifiées tout court, en posant : 


2 


2 2 2 
g ues gem -o e tes =i EE iB x 
a 


4a 2 p 
Troisième cas : z €[0,1] 


On utilise la paramétrisation suivante de l'argument z, et la fonction ¿(z) qui en découle est la 
suivante : 


Jz-z + ArcSinl/z ) 


А С@= SJ з= Сок(р) 


Ch 


_ elcos*(p))= Cos(p)Sin(p)+ ArcSin(Cos(p)) i 5їп(2р) 2 s p = зщ 


2 4 
4E). 1 [141 [Sin'(o) _1 
25219 7:5 (о) le 


D'ou l'évaluation de la fonction de Whittaker M : 


2p)-2p+x) 


WhittakerM (c, u,4k z) «2 VEE | Ge 2 Lier (z)) et J,, (Ax (z) x 


SE саа u - z) 


K 


K” 


2 1 
= WhittakerM (r, L,4K Cos? (p) x (17 + 1) cos (Sinn) 2p+ л)- ди + i) 
T УТап(р) 4 


Et l'évaluation suivante de la fonction de Whittaker W : 


С ES (Sir ele = п), (axe (2) - Соз(т(к - u), (zz) 


WhittakerW(x, ,4K 2) Wre 


Ye 
E Sg : A 
X sc) doen all z 


(тік дси 0-44): 
dE ole" usi а-и.) 


> WhittakerW (x, HA z)= \2к*е* 


K =K. 


= WhittakerW (к, HA z)= T si- AkC(z)-- de = 1) + d 


=> Whittaker W (к, HA z)= ZE саа de ES i) 


=> Whittaker W (x, HA Cos "ell T EDD 2p+ л)- дк = a 


—K 


2Kk'e T 
Whittaker Wr, u,4K Cos? (о))+ Cos et 2p)+ z) 
Tan(p) 4 


=> 


2x*e* T 
WhittakerW (к, 11,4x Cos? (о))+ Si p - Sin(2p))-- 3 
Tan(p) 4 


Estimations asymptotiques des fonctions de Whittaker donnée par NIST < Handbook of 


Mathematical functions, Whittaker functions, section 13.2.1, Uniform Asymptotic Approximations 
for large x», formules NIST 13.21.6 et 13.21.7 : 


Ces deux expressions sont également données dans le livre de 1997 (ou l'édition de 1974) 
Е W.J.Olver et W.Rheinbolt < Asymptotics and Special Functions», en chapitre 12 exercice 12.4.6 


La fonction ¿(z) se présente comme suit et sert à la construction de l'évaluation asymptotique des 
fonctions de Whittaker M et W : 


WhittakerM (=x, u,4K 2 )= Ze: É ү Lx c) z>0 
Vz + 2° mis] к? 


{= : Se 
Whittaker W (- K, 1,4x z) = É E (Z) K, (ex c) 2>0 
Б 


Une forme équivalente est la suivante : 


е Vz + z° +10842 + М2) WhittakerM (~ к, u,4k 2) J2 r2u+1) Se I (4c (z)) 


pct WhittakerW(-&, dech 2 © | e. Deckt) 
2 2 T с) ý 


Premier cas : z positif 


Je propose la paramétrisation suivante, à l'image de ce qui est réalisé par Erdelyi et Swanson, en 
1957 : 


20 Sinh(2o)+20 „y 1} 1 1 C(z) 1 Le, Sinh(2o)+ 20 
wiel g (r) MEI ce} fist =! cotanto) = e EE 


(с удо 20 


WhittakerM L K, LA Sinh? (c) x Sinh I, (k(Sinh(2o )+ 20) 


c) 20 ; 
Cosh(o) 


K 
V2 ë Sinh(20 )+ 20 „ 


Whittaker W(=-1c, ui Ak 5 (с) Sinh(o VI C T 
osh(o 


Я К, (x (Sinh(20 )+ 26) 


Deuxième cas : z très proche de 0 comme 0(i? 


Vz+z ? + Log(z+V1+z 1+ z] WhittakerM (-, u,4K z)= 2 Гои : 1) Мет, (x) 


¿(z)= D ез 
UE. us 
Ta z M5 K, (x) 


E l l+- = — WhittakerW (=x, u,4K z) = 


TA Ул os 


K 


Les formules NIST 13.21.6 et 13.21.7 sont considérées valides lorsque le paramètre x est réel 
positif. Ces formules ne semblent donc pas pouvoir s'appliquer pour le cas oü x est complexe (je 
n'ai pas trouver d'articles envisageant le cas x complexe). Toutefois appliquons sur les formules 
d'Erdelyi et Swanson, dans le seul cas oü z «0, en posant z—-y : 


Jz-zi + ArcSin(/z ) 


z e [0,1] ¿(z)= 5 = SE 


V z? -z + ArcSinhl/=z 
2 


z«0 C(z)- * 


WhittakerM (x, u,4k z) = 42 Г(2д+ ) K eg 2 (Ac (z )) |Аге(к) < S -0<z<a<l 


K 


2u 


T 
Whittaker W (- K, dr z) = (nx aM» Zeie) M 


Partons de la dernière expression de Whittaker W et prenons z=-y, il vient : 


+» + Апет) Jr +у+1ов\/у + Aal Po)-l l 
2 2 


2 


z=ë"y Iris 


| = 
WhittakerW (- c, u,Ak y) кулк Зе" |- zl, 2 (кё (y)) | 


WI 


у>0 


Or selon la formule de connexion des fonctions de Hankel et Bessel К, toujours valable quelques 
soient les paramètres, il vient : 


| = 
H » (кё (»))= Zek, " D ()= Whittaker W(- K, U,4K y) N Мак Зе e" eb) Zen e (y) 


EEN P0)-l d 
=> | = | : 
Whittaker W(- K, AK y) [ке үк, ke (y) |Аге(к) < 7] 
T Ç (y) r 


C'est exactement la formule NIST 13.21.7. Il me semble donc que cette formule puisse s'étendre 
aux valeurs complexes du paramètre x avec une contrainte sur l'argument. 


Pour ce qui est de la fonction de Whittaker M, utilisons la formule de continuation analytique 
suivante : 


WhittakerM (к, pe )= tie" WhittakerM (=x, u, z) = WhittakerM is K, dhze ) = tie ""WhittakerM (x, 11, z) 


On peut donc écrire l'expression asymptotique suivante en la reportant sur la première expression 
d'Erdelyi et Swanson : 


WhittakerM (- K, шге“) -ie ™" WhittakerM (к, u,2) x -je "^ 42 Tu) e, 2u (axe (z)) SI < 


“з Ye) 


T 
к? : 


En réalisant la méme transformation z—-y, il vient : 


(ses) чв) ny 
2 


z ery = E) 7t tmd 


PD GUY 


EE 


-o«z&a«l 


Ce qui une fois que l'on applique la formule de connexion des fonctions de Bessel J et I, toujours 
valable quelques soient les paramètres, il vient : 


Jau (axe (y))= e], (4x6(y))=> WhittakerM(-x, и,Ак у) = -ie ^" 42 r2u +1) = yh. 20) 
Y 


x ? 
SCH [Е +y 220 dr +y EE P()- 1 nl 


WhittakerM (- K, L,4K y) x42 E 


C'est exactement la formule NIST 13.21.6. Il me semble donc que cette formule puisse également 
s'étendre aux valeurs complexes du paramètre к avec une contrainte sur l'argument. 


Donc les formules NIST 13.21.6 et NIST 13.21.7 s'appliquent bien aux cas d'un argument complexe 
du paramètre x avec la contrainte de ne pas dépasser z/2. Et c'est exactement ce type de cas qui 
permet de retrouver la limite des fonctions d'onde de Bessel vers les fonctions de Bessel , en 
prenant les paramétres suivants, que l'on peut tout à fait tester numériquement : 


2 2: 2 


2 
= и=5 2= 2 О) а ева 
а 


2 = 


Construction des solutions de première espèce de l'équation de propagation de Bessel par un 
développement de Frobenius 


Ici je ne donne qu'un très bref aperçu de résultats plus complets qui seront développés dans le point 
suivant. La construction n'est envisagé que pour la fonction de première espèce. 


L'ouvrage de P.Moon et D.E.Spencer « Field Theory Handbook » réalise la construction des 
solutions de l'équation en utilisant le développement de la série autour de la singularité z=0 de 
l'équation des ondes de Bessel, et ceci par la méthode de Frobënius. On reproduit ici une partie du 
calcul réalisé par l'ouvrage. 


Recherchons la solution de l'équation autour de la singularité z=0 par un développement de la 


j=+0 
forme : Z(z)= 2" Yo. 


j-0 


Dans la nomenclature de Bócher, l'équation de propagation de Bessel prend la forme suivante : 


Z"(z)« P(z)Z'(z)+ o()z(z)- 


> Коры г20(2)= а?2* + 822? = р? 


=Т= X Az z*P(z =а?* + В?2 – р? = У Bz 
A =l A=0 i»0 
=> B E 2 = DE = 2 d = 
,=-p" B,-p' В, =а sinon В, =0 


Les coefficients se déduisent du coefficients Со de Іа série раг l'intermédiaire du système linéaire 
d'équation suivant : 


Cofo(r)=0 
C fo(r+1)+Cofi(r)=0 
C. f (r+2)+C fiir +1)+C,f,(r)= 


C fo(r+j)+C;.fi(r+j-1)+...+C,f;(r) = 


4=1 420 i»0 B,=-p° B,=B' B,-a^ sinon B,-0 
fo(r)=r*+(4-1)r B => fo(r)=(r-pPXr + p) 
f;(r)|=rA,+B;=> f(r)-0 f (r)= B? f(r)=0 fi(r)=a? f(r)=0 


Dans le cas précis de l'équation de propagation de Bessel, le système se réduit alors à la forme 
suivante : 


CGf(r)0  Gf(r«120290C,-0 

Cffr+2)}+Cf(r)=0 C,ffr+3)=0=C,=0 
Cf{r+4)+C,ffr+2)+C,f(r)=0 

C.=0 Cffr+6)+C,f(r+4)+C,f(r+2)=0 

…C,,=0 Cfifr+2j)+C,,f(r+25-2)+C,,4f(r+2j-4)=0 


L'équation indicielle donnant les racines à partir desquelles des solutions non triviales sont 
possibles est la suivante : OE (r- pYr+p)= 025n-p et һ=-р: 


Ce qui conduit à une relation de récurrence sur les coefficients pairs, les coefficients impairs du 
développement étant nuls : 


r tq f()e0 /AG()=(r-pXÁr+p)> f(r*2j) r9 2j р) 2j p) r4 2j] - p° 
fr) p" I)a? 
ee Ell. D 


net) + р? 
Er Alda ue Rl EE 


C 


= —C GC... 
P fr*2j) r«2jf-p ffe 
La solution de première espèce a donc la forme suivante : 
Ј=+0 | p 2 ^B, ,+a B,, 
B=+p => Z(z)=z"C, > B, jz? avec «(E l B, = E T В, = p SE 2 Im 
F 2) T(p+1) (p+2Ÿ =p (p+25) =p 


Avec le paramètre a=0, on retrouve le développement d'une fonction de Bessel de Première espèce : 


В? B? B? В? 
B, =l B,-- 2 В, =- 2. NB, = 2 AAA: z 
2*(p+1) 2? (p^ j) 2* j(p+j)2*(j-1)p+j-1) 
EE И раа) H _I(p+j+1) 
=С ere ( UE) IIT (p^ +1) Hens Г(р +1) 


Hon. 


I=1 


Ө e a ell 


В, 4 =L 


2) IT(p+j+1) 2) T(p+1) 
: NIC: ies (1Y E B | же (ay Ё B y s | 
xu ad DE 2 тран) 2 "Z reegen 2] ^09 


Application plus complète de la théorie de Frobënius à l'équation de propagation de Bessel 


Cas des deux racines réelles présentant un écart réel 
Revenons au cas de l'équation de propagation de Bessel гн Leila + B? -Z z(s)= 0 et 
regardons le premier où les deux racines différent d'un réel pur, il vient : 

n=p n=-p ken-n-2p ckt ei f(r)=B* Gre sinon fj(r)=0 

C (a) c (1YT(2p+1)A (p) C (n)=C, ( 1) T(1-2p)A;(- p) 


° jT(2p+j+1) jT(j41-2p) 
Го (у+17-р? 0 0 0 0 | 
B? 0 (r+2) =p? 0 0 0 
0 B? 0 (r+3)Y =p? 0 0 
А (r)= a? 0 B? 0 (r+4)Y - p 0 
0 a? 0 В? 0 >. 
x E o (r+ j—1) -p 
| 0 0 0 a? 0 В? 0 | 


Les déterminants d'indice impair sont tous nuls, comme on la déjà remarqué dans la formalisation 


par récurrence. Il reste les indices pairs et l'on peut montrer que le déterminant se simplifie comme 
Suit : 


l=j-1 1=]-1 
Comme A, (r)=0 et А, (r)=(-1Y[T(r+21+16-p?)A';(r) А, (p)=(-1) ^ GT [Gr ++ 2p) 
I=2 j I=2 j Б 5 
T IU 110520) олугоразуы) гн) 
(21+102р+1+21)= 22—55 = - 


(2j)T(2p+1+2)j 
1 


I=j- 


a, (- p)- Cy ^' C p)[ [R1+1)21+1-2p) 


BE 
1=0 1=0 


GE EE 


2; / 
— Сов Б -p)= C, d 20) a, p)= C, 2 "A A" ( p) 
[8° рор 0 0 0 0 | 
oi B? (r+4) =p? 0 0 0 
0 a? B? (r+6) =p? 0 0 
A. (r)=| 0 0 a? В? (r +8) - p? 0 
0 0 0 a? B? | 
2 (r-2j-2y-p? 
| 0 0 0 0 0 a? B? 


Il reste à définir la constante de départ Co elle est choisit de manière à retrouver, lorsque le 
paramètre a est nul, la solution classique des fonctions de Bessel. On peut montrer qu'il faut 


choisir : 
Cori) = (2) == Colra)= Bal KE 


Et dans ce cas les solutions de première espèce et de deuxième espèce, s'écrivent : 


Cyl(i=p)=(2) CI AO) C (n= 0-2). CDA C p) 


2) 29 AT (ps j+ 2, 2" jm(j«1-p) 


Los a CR s Cunt Y 


2) 4 ЛГ(р+/у+1)\2 


SE E rA A) 


2 j-0 ЛГ(7+1-р) 2 


Lorsque le paramètre a est nul, оп doit retrouver la formule du développement en série de la 
fonction de Bessel. Il suffit de considérer les déterminants : 


2 (r+2) =p? 0 0 0 0 
—Ñ B? (r+4) =p? 0 0 0 
0 0 B? (r 6f – р> 0 0 
a-0— A (r)-| 0 0 0 B (r+8f =p? 0 
0 0 0 0 B? 
m oi e (r+2j-2) 
0 0 0 0 0 0 B? 


N (=P blog" As 
zoea DEE X cy Gg 


7—0 jT(p+j+1) 2 


Ze) 1,9.) Age 5-07 (02) 7 


f j'T(j+1-p)\ 2 


C'est d'ailleurs ce qui justifie bien à rebours le choix de la constante Co 


Profitons en au passage pour exprimer les deux solutions de l'équation de propagation de Bessel 
modifiée qui permettent de généraliser la fonction de Bessel modifiée I. 


WE SE (e 2? – В? P. e= 


consistant à modifier le paramètre q° en -q*. Alors il est aisé de voir que la solution s'écrit par le 
changement de variable z->iz : 


Z (z)=J la pz) Z,(2)=J_„(a,B,iz) 
ze) 1, pas) ENK (а, (p) EE A. (p) ER 


2) 4 ЛГ(р+/у+1)\2 2) & лг(р+у+1)(2 


` Z, (z)=J_„(a, B, z)= |22) FL EE EE ap) [z] 


2 4 ЛГ(7+1- р) 2 2) < ЛГ(ў+1-р)\2 


Par choix afin de recouper avec les fonctions de Bessel modifiée de première espèce, dans le cas du 
paramètre a est nul, il est aisé de ne retenir que la valeur réelle des solutions en ôtant le facteur i? 
et i”, comme suit : 


= m zp REN А, (р) 2 MP Eo us z ; 
21,82) 5 ) H 22568) =i J (a, B,iz)=e J (a, B.iz) 
SS P j= Ад! (— z 2j ыт 
Z()=L.,(0.8,2)=| £) 2 EH =i"J_,(a,B,iz)=e ? J_,(a, PB,iz) 


Et l'on retombe bien sur la définition de la fonction de Bessel modifiée lorsque paramètre a est 
nul : 


SHARE Scheff" 


— Л\Г(ў+1+р)\ 2 
= ~ | Bz 2j-p 
Z = 0, H =] = 
-()=1,(0,8,z)=1, (B z) z 3v 2 J 


Déduisons rapidement la propriété de parité dans le cas du parametre p entier et positif pour les 
fonctions de Bessel J et I : 


J (a. B, z) = (= 1)" J (a, B, z) 


рєм -| 
I „а, В,-2)= (-1) 1, (a, B. z) 


Plan de récurrence pour la détermination des fonctions d'onde de Bessel modifiées Í 


On a vu précédemment que pour la solution de première espèce le plan de récurrence était le 
suivant : 


535 | fs. Y B? B°B,,,+a°B 
J (a, B,z)- z^C, V В, z” c. =(2) B,=1 B DN 2-2 
„(a B z) 2 €, > 232 AVEC: Gem E I(p 1) 0 2 (p+2) _ p? 2j (p+2jŸ =p 


Par extension, il est aisé de trouver le plan de récurrence de la fonction d'onde de Bessel modifiée : 


у=+о : р 1 8° В?В . —q?B | 
I (a, B,z)= z”C B,z" avec C -(2) =1 B,- B,, = = 205 
МАРО бо а E а она ере 
Avec le paramètre a=0, on retrouve le développement d'une fonction de Bessel de première espèce : 
2 2 2 2 
He В В В В, 4 =... 


=] B,=—— B, =B, „= 
i ` 2(p+l) `” £Zj(pej) ” 2 EVENT EH CET ER 


B” 1 (BV Dos Fie, n IP+j+1) 
u | -(2) NT (p+ +1) RAGE Г(р+1) 


saf) Leur (2) 
# (2) AT(p+j+1) ° (2) Г(р+1) 


I DLE) S 1 |22) Е 5 1 (2) ris 


2, Z (pe j«l) 


Cas des deux racines réelles ou entières identiques : le seul cas ici se produit à p=0 


Les deux racines de l'équation indicielle coïncident, la première solution reste identique : 


сй! c0) 57 010) ae (a, nal SONO) 


722 яг(/+1) Z. (y 


B 4 0 0 O 0 | 
a BP 16 0 0 0 
0 aÏ 8? 6 O 0 
A',(0)=|0 0 a B? g 0 
0 0 0 oi B . 

e (2j-2Y 
lo 0 0 0 - 0 D 


20) testa) (=o 06а) | эз j 


r =0= El Los(s) fa. ëch 5 US E ШҮ» ] 


= (jy 


Puisque 20=0, p=0, et que les déterminants d'indice impair sont tous nuls, il vient : 


Ikon 


Z, (z)= Log(z)J (a, В,2)+ > Zu (GI dÉ, 


p=0>3 A,,(r)=(-1YT ee V A',(r)=(-1) ПУХ A' (r) 


ve 

GN 
ll 

= 


si N 2^0). Cay Rea ara Есе | 


Posons g(r,l)=(r+1IY M}, (r) cofacteurs ligne i colonne k de A (r) 


1 


aen E) era 8609) ру, 0) ы Ai 


r 


8er.) з). 98,0) ш) L (r)+(r +4)м/ (r)+...+(r+2(j-1))M L, (r)! 


Or ôr "m 
= SEXO = pu (sae) nt DM C] e ti) ent) G7 01,0) 
pe stierfs De loco BOUET, Si: 
юн Dn [Oo] 

eus ku RDO TI" 
=> Z,(z)= Log(z)Z (z)+ 2 o pO 
-2(-1} 2 (jn? ^' (0) - 
=> Z,(z)= Log(z)Z ( D y 1) О p 20 2w (0), a 2A', DI 
= Z,(z)= Log( EXC iy Em to) резу. 


äm i=j Ј=+о 


De plus У - Y - 1$ 7 (02 лові) É C = УК дна) 


1=1 m 1=1 m jA =| 


= 


On peut s'arrêter à cette expression, mais il est d'usage de recoler avec la définition sous forme de 
série de la fonction de Bessel de deuxième espèce Y0 (lorsque a=0). Et dans ce cas il d'usage de 
rajouter un terme à la seconde solution, comme suit : 


Y (a, B, z)= 202, СӨ) TE) 


nlepa)? nl Zille Слу 27, LS gea 0 


1=1 1=1 


On donne donc l'expression formelle pour la construction de la seconde solution : 


220070 pa RAAN] 


y constante d'Euler у = 0.577216 


=+00 т 2т m- m 
Y. (a. В, z)- КЕ 3) (a.f, z)- bx =) B PO LA, D» 
1 (m!) 2 I ml 
|^ 4 0 0 0 0 | 
aÇ BP 16 0 0 0 
0 a р & o 0 
A',(0}=1 A',(0)=Del 0 0 oi pos 0 «-m-]l termes 
0 0 0 aÏ B ET 
" (2m - 2) 
[0 0 0 0 0 uw P 


^' (0) = Dep’ ]- 8? x. (0)= Del É, ML 


а 


M 70) le mineur du déterminant A, (0) 


Le mineur Mu, (0) est le mineur du déterminant A',(0) obtenu en ôtant la ligne let la colonne 1+1. 


Une conjecture que l'on peut formuler se rapproche de la formule de construction des fonctions de 
Bessel traditionnelle : 


"pens Ale, B. z)Cos(vz)- J (о, B, z) 
(2,8,2) та Sin(vz) 


On montre que la formule est bien valable, mais il faut auparavant avoir établi deux choses : 
- la formule de connexion des fonctions Y tilde (construite à l'aide des fonctions de Whittaker) 
- la formule du décalage entre les fonctions Y et Y tilde. 


Dans le cas où le paramètre a est nul, le déterminant ayant la forme, tous les mineurs calculés sont 
automatiquement nuls, et il vient : 


= Z (z)=(8 z)= > T: 


(= = 2. y+ Log LE | 2)- E - Bl D ” (0)+A'„ (0) Y J 


l=1 


Ba 
KE y constante d'Euler 


B 4 о о 0 0 | 
0 д 16 0 0 0 
оор 6 O 0 
A'(0O)=De 0 о о P 8 0 =?" MP" (0) le mineur de A'„(0) 
оо о 0 B T 
" (2m - 2/ 
10.0 0000 f | 
[0 16 0 0 0 | ва о 0 0 | 
ов eco 0 0060 0 0 
= M"(0)- Del 0 y P id ~ Q M},(0)= Det G ust P = 
' оо о g 0 | 0 0 0 pB 
" e (2m-2Y = . (2m-2) 
оо 0 0 B | ооо 0 B A 
B 4 о 0 0 | 
0 f? 16 0 0 
M” A [0)= Det : : g P ' =0= Vlefl.,m-1}M},,(0)=0 
0 0 0 O 0| 


Ce qui justifie la transformation opérée sur la seconde de solution dans le cas du paramètre a est 
nul pour retrouver l'expression usuelle de la fonction de Bessel de deuxième espèce Y. 


Cas des deux racines entières distinctes ou de deux racines réelles présentant un écart entier 


Là encore je reproduis une partie du résultat de la première et seconde solution tel qu'exprimée 
dans l'ouvrage de Moon et Spencer. Cela permet également de s'assurer qu'il n'y a pas d'erreur de 


calcul. 


La première solution reste identique : 


p» (2 ECC (EJ DAO zesg (2) S CI 8 Oz f 


2) n! 2) XY = Л 
8° (n+2) =n? 0 0 0 0 
a? B? (n+4) =n? 0 0 0 
0 a? B? (n+6) =n? 0 0 
A. (n)= Den 0 0 a? p (иву =n? 0 
0 0 0 a? p Е 
Та ii (n+2j-2)Y =n? 
0 0 0 0 0 a? B? | 


Pour la seconde solution partons de la formule générale de la théorie de Frobénius et appliquons 
là au cas de l'équation de propagation de Bessel : 


z=0 љ=-р A, p)=0 k-2p Tous les termes impaires des sommations s'annulent !! 
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Qj) WW 
Pour j»p A,(=p)=(=1Y*” "OH 


De plus Pour j=p 


Pour établir ces trois dernières expressions des déterminants complets vers les déterminants réduits 
on effectue un calcul similaire à celui du cas des racines réelles distinctes, à ceci près que 
maintenant les racines sont entières et que le produit désormais exclue le terme qui annulerait le 
dénominateur, comme suit : 
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Démontrons également une formule sur les déterminants réduits, dont on remarque qu'elle est 
incorrectement traduit dans l'ouvrage de Moon et Spencer, tout en n'occasionnant par d'erreur de 
calcul ultérieur car elle est correctement utilisée. 


La formule sur les déterminants complets est donné ici telle quelle lorsque l'écart des racines 
indicielles est entier : 
Asa) 9 Als Jä (n) А, as p)7 ^s, C p)A,,(p) 


Substitution toutes les expressions trouvés précédemment pour établir la méme formule sur les 
déterminants réduits : 
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L'expression exacte sur les déterminants réduits dans le cas du développement de Frobenius pour 
la fonction d'onde de Bessel est donc : 


A' (C p) А", (C p)A' (p) 


Pour ce qui est des dérivées des déterminants, il vient : 
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La substitution de toutes ces sous expressions dans l'expression principale de la seconde solution 
donne : 
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Les deux dernières étapes du calcul consiste à décaler la deuxième sommation sur j pour la faire 
partir à 0 et de même que dans le cas p=0, la seconde solution seras légèrement modifiée pour 
recaler avec la définition classique des fonctions de Bessel de deuxième espèce Y d'ordre entier. 
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En introduisant ensuite la fonction d'onde de Bessel de deuxième d'ordre entier Y, comme suit : 
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On donne donc l'expression formelle pour la construction de la seconde solution (p->n,j->m) : 


П j=+œ —1)/А'. 2j 
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B^ (n+2f =n? 0 0 0 0 
a? B? (п+4) =n? 0 0 0 
0 a? B? (n+6) =n? 0 0 
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[8° (2-nY - n° 0 0 0 0 | 
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Le mineur Min, (-n) est le mineur du déterminant А'„+„(-п) obtenu en ôtant la ligne 1 et la 
colonne Ї+1. 


Les deux déterminants introduits sont de taille (jj) et peuvent se construire comme fonctions des 
indices j, n et p ainsi que des deux paramètres a et B. On peut tous deux les caractériser par une 


récurrence et une valeur de départ : 


[8° (2+p) -n° 0 0 0 0 
a? B? (4+ pf =n? 0 0 0 
0 a? B? (6+ p) =n? 0 0 
A' (a, B,n, p)- Det 0 0 a? p? (8+ p) -n 0 
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Il est aisé de montrer qu'il peut se calculer par une loi de récurrence : 


7 (a, B,n, p)= D'A, (a, p,n, p -2)- a^ (2+ p+n)2+ р-п)л',, (a, B,n, p + 4) 
A, (a, B,n, p)- 0 A, (a, В,п, р)=1 


Plus précisément pour les deux déterminants suivants, оп а la récurrence et son point de départ : 


N (a, В,п,+п)= B? N', (a, B,n,2+n)-a?*(2+n+n)(2+4p-n)A"'„_,(aœ, B,n,4 +n) 
A', (a, В,п,+п)= BON. (a, B,n,2 + n)- 4a’ (1 + n)A' , (a, B,n,4+ n) 
A, (a, B,n,+n) = 0 A, (a, В,п,+п) =] 


Ce qui permet de calculer par exemple les quelques premiers déterminants non triviaux : 


А (a, B,1,+1)=q?" => AN (a, B,2,+1)=q? 
A', (a, 8,2,2)= B* -12a? 


л", (a, B,2,#2)=q?*A"' (a, В,2,2+ p)- 4a? (1£ 2)^', (a, B,2,4 € p)= B^ - 4a? (152) " 
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La formule redonne le développement en série des fonctions de Bessel de deuxième espèce Y : 
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La méme conjecture sur la seconde solution n'est ici vérifiée numériquement que pour la valeur 
р=1: 


m J (a, B. z)Cos(vz)- J (о, B. z) 
Y (a. B. z) M Sin(vz) 


On montre que la formule est bien valable, mais il faut auparavant avoir établi deux choses : 
- la formule de connexion des fonctions Y tilde (construite à l'aide des fonctions de Whittaker) 
- la formule du décalage entre les fonctions Y et Y tilde. 


Résumons donc les expressions obtenues des développements de Frobenius, pour les premières et 
secondes solutions, de l'équation d'onde de Bessel : 
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Et pour le cas a=0, redonnant les classiques fonctions de Bessel, J et Y : 
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Quelques expressions concrètes de développements limités autour de l'origine 


Voici le développement limité à l'ordre 8 des fonctions de propagation de Bessel d'ordre 0 de 
première et deuxième espèce : 
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y constante d'Euler y = 0.577216 


Voici le développement limité à l'ordre 10 des fonctions de propagation de Bessel d'ordre 1 de 
première et deuxième espèce : 
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y constante d'Euler y = 0.577216 


Voici les développements limités à l'ordre 10 des fonctions de propagation de Bessel d'ordre p et -p 
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Voici les développements limités à l'ordre 10 des fonctions de propagation de Bessel d'ordre 1⁄2 et 


-1⁄2 : 
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Il est à noter qu'il semble y avoir une erreur dans le report des formules dans l'ouvrage de Moon- 
Spencer « Field Theory Handbook » à la page 189. 


Plan de récurrence des solutions de Frobenius régulières et non régulières de l'équation d'onde 
de Bessel autour de z=0 pour des valeurs entières de p 


Le plan de récurrence de la solution régulière est le suivant : 


j=+œ p =0 B,-1 
J =z"C В, z?! C == ЕЕ д\ | 
j (a, S z) < Me p) 2 Gi TH dr p) 2T(p + 1) Ps БЕ 1p EJE US, F q” B + а B, j, =0 


Par exemple pour p=0, il vient : 
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Le plan de récurrence de la solution non régulière est établi selon la forme choisi du 
développement autour de z=0 auquel on applique l'opérateur différentiel correspondant à 
l'équation d'onde de Bessel : 


th Late latz + B? Eat 0 


Prenons pour le cas p=0, une solution de la forme suivante et appliquons lui l'opération 
differentielle de l'équation d'onde : 
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Dans ce cas р=0, la récurrence commence au coefficient A» tandis que la récurrence pour p» 0 
commence à Ao. 


Pour le cas p>0, prenons une solution de la forme suivante et appliquons lui également l'opération 
differentielle de l'équation d'onde : 
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On note que la récurrence précédente ne permet pas de calculer le coefficient Ao (où j=-1) et qu'il 
nous faut le déterminer par son expression littérale. 


La constante À, est un déterminant calculé comme suit : 
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Pour ce déterminant, il est aisée de montrer qu'il peut se calculer par une loi de récurrence : 
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Ou p est un entier. Plus précisément pour les deux déterminants suivants, on a la récurrence et son 
point de départ : 
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Ce qui permet de calculer par exemple les deux premiers déterminant : 
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Autrement écrit, on a également : 
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Par ailleurs comme nous l'avons remarqué la récurrence ne fonctionne pas concernant le terme A , 
pour cela on peut revenir à l'expression littérale que l'on a calculé précédemment : 
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Cette expression littérale donne à l'aide de Mathematica, les expressions suivantes : 
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On peut donc maintenant énoncer la récurrence complète des termes du développement de la 
seconde solution de l'équation d'onde de Bessel : 
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Il est clair que cette récurrence ne rend compte dans l'expression choisie : 
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Le premier terme qui est un terme de décalage de la solution régulière suit la récurrence de cette 
dernière, indépendamment des deux autres termes. 


Pour la fonction d'onde de Bessel de deuxième espèce construite par Moon et Spencer qui s'écrit 
comme suit : 
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On a le plan de récurrence suivant : 
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Il se dégage du plan de récurrence, une valeur limite intéressante à calculer : 
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Nous allons maintenant essayer de déterminer le < décalage > en les solutions de deuxiémes espèce 
Y proposées par Moon et Spencer et les fonctions Y construites avec les fonctions de Whittaker, 
notamment à l'aide des comportements asymptotiques de ces fonctions. 


Comportement asymptotique des fonctions d'onde de Coulomb pour n proche de zéro avec un 
paramètre p entier 


Pour n proche de zéro, le comportement asymptotique des fonctions d'onde de Coulomb est le 
suivant 
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Cela correspond au comportement asymptotique des fonctions d'onde de Bessel telles que calculées 
d'après les fonctions d'onde de Coulomb : 
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Il convient de préciser que le comportement asymptotique des fonctions d'onde de Coulomb et donc 


2 
a>>0,B tq P ага 
de Bessel est atteint dés lors que le rapport des paramétres а est très petit. 


Il est possible de valider le comportement asymptotique de la fonction d'onde de Bessel de première 
espèce d'après la formule littérale du développement de Frobënius : 
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En effet le calcul direct des déterminants est possible dans le cas a<>0 et B=0 (ou proche de zéro 
par rapport à a), qui correspond au comportement asymptotique recherché. Dans ce cas les valeurs 
des déterminants sont les suivantes : 
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Injectons l'expression du déterminant dans le développement : 
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On retrouve exactement le comportement asymptotique recherché pour la fonction de première 
espèce. 


Regardons maintenant le comportement asymptotique de la fonction de deuxième espèce. Calculons 
la forme de la fonction de Bessel asymptote et commençons par le cas le plus simple lorsque p=0 : 
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En partant du développement de Frobénius de la fonction de deuxiéme espéce, et en posant 
formellement B=0, il vient : 
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Les déterminants se calculent comme précédemment en posant p=0 : 
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Les mineurs ou cofacteurs ne sont pas tous nuls en revanche, ils s'annulent pour l'indice m impair, 
on a donc : 
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Par un calcul effectué à l'aide de Mathematica, on peut vérifier l'égalité suivante sur le terme 
polynomial du développement de Frobénius, par identification des termes limités dans la série : 
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On reconnatt le développement en série de la deuxième partie polynomiale de la fonction de 
Neumann Y. Il vient alors : 
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Lorsque n est un entier non nul, le développement de Frobenius donne : 
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Les déterminants se calculent comme suit : 
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Par ailleurs dans l'expression des fonctions de Bessel de deuxième espèce (Neumann), il faut 


distinguer les demi- ordres selon que l'ordre est paire ou impaire : 
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Lorsque l'ordre est impair on voit tout de suite que le terme de première espèce est nul. D'autre part 
on peut montrer que l'addition du terme de puissance inverse et du terme polynomial est 
proportionnelle au développement de la fonction de Bessel de première espèce, d'argument ax°/2 et 
d'ordre -p/2, plus précisément du terme proportionnel déjà évoqué dans l'étude du comportement 


asymptotique des fonctions de deuxième espèce : 
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calcul à l'aide de Mathematica 


2 2 p+1 "x | á 2 
~ Ц x 
dëch да) (5) TZ |=) 


Donc on conclut qu'il n'y a pas de décalage pour les ordres impaires. 
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Il reste le cas des ordres paires. La fonction d'onde de Bessel s'écrit comme suit : 
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D'autre part la fonction asymptote a la forme suivante : 
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On montre à l'aide de Mathematica, que le terme de puissance inverse et le terme polynomial dans 
le développement est proportionnel : 
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Unifions les formules avec le cas p=0, il vient : 


TT з PORTE) 


1 1 23 1 1 176 1 1 1 
A(0)=0 All)=- А(2)=—=1+— 3) = =1 A(4)= =1 
(0) () 1 e) 1 t3 6) 1x3x5 *3 5 (4) 1х3х5х7 73 39 
aee ета PN ЯЕ A7 = "m 
315 3 5 9 3465 3 5 9 11 45045 3 5 9 11 13 
On émet donc l'hypothèse : av 1 
o 21+1 


Voici quelques graphes de fonctions pour des ordres pairs, qui illustrent parfaitement la valeur du 
décalage à opérer sur les fonctions asymptotes étudiées : 


a = 70.5, B = 1.5, p = 2, nSeries = 20 
be 5) Веѕѕе Маме Series S2Function(70.5',1.5',2, x ,20) 


22 „70.522 d A) 


“х2 
= le x ) 
2 2 


"x? 
"m ya (== |+ 
2 2 


0.50 


.2 | 
еза) 
: 2 


л 


2 52 
"und E Ja 
2 


| 


70.5 


2 


T 
= vy [= x ) 
2 т 


a = 70.5, B = 1.5, р = 4, nSeries = 20 
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Dans ces trois graphes le rapport B°/a n'est pas assez petit: 


0.4 


а = 70.5, B = 3.5, p = 2, nSeries = 20 


pe 24 ) BesselWaveSeries SZFunction(70.5,3.5',2,x,20) 


22 70.522 d 


den E | 
ME CU =< ss 


E kel 


Si A (222 ) 
2 2 


06 


a = 70.5, B = 3.5, p = 4, nSeries = 20 
оь 


bei 2 (nx BesselWaveSeries S2Function(70.5'.3.5'.4,x 20 


ez| gas?) 
на) 
=, ВЕ) 
– ab 
» La EI 


edi gl 


a = 70.5, B = 3.5, p = 6, nSeries = 20 
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Dans ces quatre cas pour p impair prenons le rapport B°/a suffisamment petit : 


a = 70.5, B = 1.3, p = 5, nSeries = 15 
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Dans les trois exemples suivante le rapport p°/a est suffisamment petit : 
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a = 70.5, B = 1.2, p = 12, nSeries = 20 
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Conclusion : on obtient la forme asymptotique des fonctions d'onde de Bessel de première et 
deuxième espèce définie par le développement de Frobenius et la forme Y tilde: 
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Hypothèse de construction de la seconde solution 


Au vu des définitions des fonctions d'onde de 
asymptotique, il vient : 


Bessel suivantes, et de leur comportement 
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Comme on a mis en évidence un décalage entre les deux solutions de deuxième espèce lors de 
l'étude du comportement asymptotique pour a très grand, il vient : 
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Commençons par Іа valeur р=0, et retranchons le terme divergent еп logarithme, dont nous 
connaissons la valeur, des fonctions de deuxième espèce, il vient : 


Yoa, B,x)=Z Logs) У" lapa) Die, B.x)= Log(x) +V" (a. 8x) 


or Lin leu s) 2 Log(s)) = {у + Lou Ë) et LimJ (z, B,x)=1 

De plus B.x)- Zil - Ze À EE Я 

sou a, ede Rie 2н) 20 bd (A oet 
-F (a, -H DCE 13) LE 


=> Hypothèse Y,(a. B, x) = Y (a, В,х)+ о | _ DCE —i Ch (o B, x) 
T 


4a 2 4a 


Cette hypothèse est cohérente avec le comportement limite a->0 : 


x20 x 


2 
x= uL PTS pig 3 d Log 29 |+Rd y| 277 le 2 (3p*-1) 
B 2 x>0 (у B 2 4a 3p 


= a 0 (6) (25 га) = (3p? -1) 


DEER SÉ Jj = en +6и?)= (1+. виби -1)) 


Cette hypothèse est également cohérente avec le comportement limite a très grand puisque : 


Y (a. p, x) = e poe Hro E (1) nta 
T 4a 2 


Et di =_-y-2Log(2)= Y,(@, B, x) « (о, B, x) EZ +y + tos) fefe 
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G 
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Nous avons trouvé le décalage exact en fonction d'onde de Bessel de premiére espéce, pour les deux 
seconde solutions définies l'une par la série de Frobénius et l'autre par les fonctions de Whittaker 
ou de Coulomb, à savoir : 


Y o. B.x)- Èa. pr All Z) DCE Elon, B.x) 


Le résultat est également confirmé par la visualisation des graphes suivant, mettant en évidence la 


formule du décalage : 


son 


#1; 
= 1; 
-e; 
— BesselWaveSeriesS2Function(a, В, p, x, nSeries) 
— BesselYWaveFromWhittaken(a. B, p, x) 


[ees jefo 3 ))pcoseuwaverrommataxenta sa D 
— BesselYWaveFromWhittakerW(a, 8, p, x)+ 


nSeries = 15; 
0 


—— BesselWaveSeriesS2Function(a, B, p, x, nSeries) 
— BesselYWaveFromWhittakerW(a, B, p, x) 


EI E) rol ue )|)=esceuwaverromninakeruca s.p. x) 
— BesselYWaveFromWhittakenW(a, 8, р, х) + 


n 


Par le même passage à la limite sur les ordres entiers positifs, partons du développement limité 
suivant : 


a(r + koe Pe tont (=n)J„(@. B. x) — E 19 E 
rp" 2 x , 


p" m-0 m! 2 
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T 
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Par nature le terme logarithmique est proportionnel à la fonction d'onde de Bessel de première 
espèce. Pour le terme inverse, il s'agit d'un polynôme de degré négatif x” et de degré positif x"?. Et 
nous allons supposer par hypothèse la forme suivante du décalage : 


ев) ëch A r+ Log] lais pell (n), (o. B.) 


л B?" 
Logla?)+ ordy" -iE aan 
g(a,B,n)= y + 2 S 


Log(a?)+ LUE 200 J- Ala, B.n) 


4a 
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4 


Posons également Blo, В,п)= А, (=n) Аа. Вл) = Л! (-п)С(о, В,п) 


À l'aide de Mathematica, nous avons calculé la limite de la fonction Y tilde, il vient pour les ordres 
n=1,2,3 : 


Logla?)+ indui ab 


| A(a,B1)=0 Ale, B,1)= B? 


п=1 glo B,)= y ^ ` 2 
Logla?)+ 2 DCE -i £) - A(a, 8,2) 
n-2— g(a,B,2)=7y+ 7 2 
A(1,1,2)= ` A(1,2,2)= = A(1,3,2)= М A(2,1,2)= = 4(2,2,2)=1 д(23,2)= = 
Or A',(-n)=A(a,B,n) A(a,B,2)=4a° + В“ 
8a? Aud x2a? 
A(1,1,2)=5 A(1,2,2)=20 А(13,2)= 85 ai Ale, B,2)= SC = EIERE 
A(2,1,2)=17 A(22,2)=32 A(23.2)=97 
> Bla, B,2)= 4a 
Logla?)+ 2 SUE -i £) - Ala, 8,3) 
n=3> g(a.B)- y * " ` 
; _ 240? f) _ V x3a? B? 
A(a, B,3)= mE Ala, B.3)- B^ (16a? + 8*)= Aa. P A(a,B3) л(о, B.3) 


> Bla, B,3)=12a? B? 


Soit les trois formes suivantes du décalage que nous allons tenter d'étudier sur les ordres entiers 
supérieures : 
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2 4a 
ou bien 
Eland leh: SE Nc A pan) Lod E )- veh ua 
ou bien 


nia Pla pay =L CODE va DCE HIE 
л В 4а 2 4а 


I-p-l I=p 


De plus A, (-2p)= Пого в") As, CQ p) в Па? + 8°) 


1-0 I=1 


Pour n=], il vient : 
Hate: ese iro E neo Jas E 
л 4a 4a 


Il n'y a pas de décalage dans le cas a->0, puisque : 


4 1 2 2a? 2 
a 0 {зе nd: £) = (n? -1) 


Et dans le cas a trés grand, le décalage disparaît également puisque les fonctions de deuxième 
espèce d'ordre impaire augmente comme des puissances positives de а supplantant le terme 
logarithmique. 
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Pour n=2, il vient : 


Yn. x) (а.з) a ZÉ ) Ala, 8.2) и 3-2) Ja. (2). (px) 


4a 2 2 


= 4a + B° i 4a? 3 ng 
Y,(a. В,х)= Y. (a, B, ch DE). IP (2-7. ,x) 
Dans le cas a->0, le décalage s'annule, puisque : 
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T B 2 4a 
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= Y,(a,B,x)= Y (a. B.x) 


Et dans le cas a trés grand, le décalage s'écrit : 
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л p 4a 2 


v(2) -24(2)- 2Log(2)- y =2-2L0g(2)-7 = Y, (a, B.x) « Elo, В,х)+ Lo É ОЕ 
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Ce qui est encore le décalage trouvé dans le comportement asymptotique a trés grand. 


Pour n=3, nous avons : 


: СЙ 16 2 p* 2 12 2 p? 
A',(-3)= B*(16a2+ 8*)= Y. (o. В, x)= % (а, В, х)- = ta a Re | (22% 


a> л a 


a0= йй dl - (0, 8, x)+ Lim 4. SEN Re d 2. Jy (а, В,х)= (а, В, x) 


Га valeur du décalage est confirmé раг les graphes suivants : 


nSeries = 15; 
— BesselWaveSeriesS2Function(a, 8, p, x, nSeries) 


KE £4 Vo, i UE DI 
— BesselYWaveFromWhittakerW(a, B, p, x) + 


т 


— BesselYWaveFromWhittakerW(a, B, p, x) 


p=2; 
nSeries = 15; 
05 1.5 


— BesselWaveSeriesS2Function(a, B, p, x, nSeries) 


erts (mne Вр.х\ 


—— BesselYWaveFromWhittakerW(a, B, p, x) + 


ngt 
— BesselYWaveFromWhittakerW(a, B, p, x) 


p-3; 
nSeries - 20; 


— BesselWaveSeriesS2Function(a, B, p, x, nSeries) 


(16 a2+8*) Е! Ep Ru UE x) 


—— BesselYWaveFromWhittakerW(a, B, p, x) + 


ngt 
—— BesselYWaveFromWhittakerW(a, 8, р, x) 


Étudions le décalage pour les ordres 4, 5 et 6 : 
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B(11,4)=216 B(12,4)=576 4(1,3,4)= 2136 
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Pour n=4, il vient : 
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g (ка + в* зва? + в*){, (B2), лава? +в") sg 
Y, (a, 8,х)= Y. (o, B, x) л В Гов 4a Ц We" 8^ 3602 + p^) Ке ү 2 "as J (e, B, x) 
_ 2 2 
po) sos) seb ar) moet EE SUE EI Jala, B, x) 


Pour n=5, l'expression du décalage entre les deux fonctions d'onde de Bessel de deuxième espèce 
donne : 


Kloe, B. x)-- 
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2 p2 E 2 2, oi 2, oi B? p? 
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Pour n=6, nous avons : 
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Ys Bos) sfr) stat vci et ose? mets P) QUE Ch J (e, B.x) 
T 


Les graphes suivant confirment la valeur des décalages pour n=4, le décrochage entre les deux 


courbes est lié au nombre de termes restreint pris dans la série dont la convergence s'avère 
beaucoup moins rapide avec les valeurs croissantes de l'argument : 


B =1; 
p=4; 
nSeries = 15; 


, — BesselWaveSeriesS2Function(a, B, p, x, nSeries) 


24 a? (8 a2+8*) (L5 iB? 
-Ref 9) (.18* 4 ) 
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—— BesselYWaveFromWhittakerW(a, B, p, x) + - 
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—— BesselYWaveFromWhittakerW(a, 8, p, x) 


a=5.5; 
2 В=1.3; 
p-4j 
nSeries - 25; 


— BesselWaveSeriesS2Function(a, B, p, x, nSeries) 
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бе Ves 2i 
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— BesselYWaveFromWhittakerW(a, B, p, x) 


Pour n=5 : 


@ = 6.5; 

B = 2.3; 

p =5; 
nSeries = 15; 


— BesselWaveSeriesS2Function(a, 8, p, x, nSeries) 


40 a? (28 a2+84}416 a2+8*) (64 42:64) [ [3 490) É BesselJWaveFromWhittskerM(a,B.p,x) 
— BesselYWaveFromWhittakerW(a, 8, р, x) + " 
— BesselYWaveFromWhittakerW(a, B, р, x) 


a = 6.5; 
B = 2.3; 
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nSeries = 25; 


— BesselWaveSeriesS2Function(a, 8, p, x, nSeries) 
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Ha 


— BesselYWaveFromWhittakerW(a, 8, p, x) 


150 
100 
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B=2.3; 
" £ 1 D p = 6; 
1.4 1.8 2.0 22  nSeries = 15; 
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-100 
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a = 6.5; 
822.3; 
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—— BesselWaveSeriesS2Function(a, B, p, x, nSeries) 
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=p" 


— BesselYWaveFromWhittakerW(a, B. p. х) + 


— BesselYWaveFromWhittakerW(a, B, р, x) 


Revérifions pour le cas limite a très grand que l'on retrouve pour les ordres pairs, le décalage déjà 
établi, au moins pour les ordres n=2 (déjà réalisé), n=4 et n=6 : 
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À 2p 2 2 
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Ce qui est le décalage déjà trouvé dans le cas a très grand pour les ordres pairs. Pour les ordres 
impairs il n'y a pas de décalage puisque : Napa С@р+1ў= Ge 


Lim Ala, B,n)= 0 Я ; Ka 
Pour le cas a->0, comme «^ , le décalage s'annule comme déjà vu auparavant en ce 


qui concerne le terme logarithmique. 


Il reste maintenant à trouver la formule exacte du décalage pour toutes les valeurs de l'ordre n. 


Visser da - [ Geer ЕЕЕ 


1=р-1 l=p 

avec A^, (-2p)= (21+1) a? + g*) A'apai (-(2p+1))= В? (16/2a2+ 8*) 
1=0 l=1 

B(a, B,0)= 0 

Bla, B,1)= 0 

Bla, B,2)= 4a? 
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Bla, B,5)= 40a? B? so? + В“) 
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Bla, B,7)= вда? p? да? + 8* (оба? + p*) 

B(a,B.8)=112o 2(42240a +8464a*B* + 200a? B + в") 

(о, 8,9) = 1440? B ^ (ssa? +8 * Jas84a^ 12202? p* + B 8) 

Bla, B,10)=180a 2( 90800644 8 +1971712a*B* + 562240 485 +448а2 p? + В!) 


Co 


La première constatation concerne le premier terme multiplicatif : 


Ba, B,n)=2n(n-1)C(a, В.п) 
С(а, B,0)=0 

Cla, B1)=0 

Cla, B,2)=a? 

Cla, B,3) =a? B? 

Cla, B,4)=a? 8a? +85) 
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Cla, B.10)- o s +1971712a*B*+56224a* B$ +44802B1 + p) 


Pour trouver la forme exacte du décalage, il faut utiliser une formule explicite du développement 
limité autour de 0 de la fonction de Whittaker W. Cela nous permettra de pouvoir établir la valeur 


de la limite suivante : 


Y (a, Jer toss x)A' , (a. B, -p)xJ,( a, ,x)- ZA SR 


Lim 
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Commençons par l'ordre 0: Lim Y (a, B,x)-2 Log(x)xJ (o; B,x)- D'aprés les développements 
x! Л 


limités de Іа fonction de Whittaker W et M, au moins à l'ordre 0 : 
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Pour n=0, la limite s'écrit : 

е 2 T 2 1 . B? 

A (0)=1— Y,(a B.x)—-— Log(x)= Y (0, B,x)- — Log(x)« = Іов(а)+ Re y| —+i— ||+2y 
T T T 2 4a 


Dans le même temps la limite avec la fonction non tildé est la suivante 


EEN 


D'ou le décalage : 
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Selon des formules du développement limité autour de 0 donnée par le site de Mathematica tant de 
la fonction M que W (formules 07.44.06.0013.01 et 07.45.06.0017.01), on a : 
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Ce qui peut s ‘écrire selon trois termes de développement (inverse, logarithmique et polynomial) : 
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Le terme inverse est une expression polynomiale de degré négatif -u 1/2 et positif u-3/2. Le terme 
logarithmique est directement proportionnel à la fonction de Whittaker de première espèce M. On 
va retrouver cette méme structure dans la fonction d'onde de Bessel Y tilde. La limite qui devient 
intéressante à calculer est alors la suivante : 
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L'application du développement limité de la fonction de Whittaker W à la fonction d'onde de Bessel 
Y tilde : 
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Conduit également à la séparation en trois termes : 
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Donnons juste le terme logarithmique plus ou moins facile à calculer : 
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Au passage on confirme facilement le résultat suivant sur l'expression des fonctions Gamma : 
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soustraction du terme logarithmique et inverse. On peut donc continuer l'évaluation de la limite : 
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Le résultat particulièrement important qui va nous permettre de calculer le décalage exact entre les 
fonctions de Bessel Y est donc le suivant : 
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Comme Іа limite similaire ауес Іа fonction de Bessel Y de Moon et Spencer s'écrit : 
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Continuons un peu ce calcul : 
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Vérifions pour quelques valeurs de p que cela confirme les résultats déjà obtenus. Pour p=1 c'est 
bien le cas : 
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On retrouve précisément tous les résultats énoncés précédemment. 


Une fois déterminé le décalage continuons d'étudier les développements limités des fonctions 
d'onde de Bessel à partir de ceux des fonctions de Whittaker. 


Le développement limité de la fonction de Whittaker M s'écrit également 
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L'application de la première formule pour Whittaker M pour la fonction d'onde de Bessel J donne 
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Il convient de remarquer que cette expression est bien réelle 


Reprenons l'expression du développement limité de la fonction de Whittaker W. Le terme inverse 
pour la fonction de Whittaker W se calcule comme suit : 
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Ce qui donne pour le terme inverse de la fonction d'onde de Bessel Y : 
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Et finalement le terme polynomial du développement de la fonction de Whittaker W donne : 
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Donne pour la partie polynomiale de la fonction d'onde Y tilde : 
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En résultat les deux développements des fonctions d'onde de Bessel J et Y tilde : 
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Le décalage sur les fonctions d'onde de Bessel Y et Y tilde étant celui-ci 
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Оп еп déduit tout de suite le développement limité explicite de la fonction d'onde de Bessel 
construite par Moon et Spencer : 
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Voyons maintenant comment se présente les développements limités des fonctions d'onde de Bessel I 
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En utilisant la formule de connexion pour la fonction d'onde de Bessel modifiée I, il vient 
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Passons maintenant au développement limité des fonctions de Hankel telles que définies : 
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Comme оп а les trois termes du développement de Іа fonction de Whittaker W : 
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En les appliquant sur la formule de la fonction d'onde de Hankel H2 : 
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Soit pour les deux fonctions d'onde de Hankel : 
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En appliquant les deux formules de connexion pour les fonctions d'onde de Bessel I et K, il vient : 
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Formules de connexion des fonctions d'onde de Bessel Y construites par Moon et Spencer 


Je complète ici les deux conjectures sur les formules de connexion à l'ordre O et 1 pour ces 
fonctions, une fois que nous avons pu déterminer le décalage et les fonctions d'onde de Bessel Y et 


Y tilde. 
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Avec la formule de décalage entre les fonctions de Bessel Y : 
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